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ABSTRACT 

We  present  a  collection  of  algorithms,  all  running  in  0{n^  log  n  log'n)  time, 
for  constructing  a  skeleton  representation  of  a  suitably  generalized  'Voronoi 
diagram'  for  a  ladder  moving  in  a  2-dimensional  space  bounded  by  polygonal 
barriers  consisting  of  n  line  segments.  This  diagram,  which  is  a  2-dimensional 
subcomplex  of  the  3-dimensional  configuration  space  of  the  ladder,  is  introduced 
and  analyzed  in  a  companion  paper  by  the  present  authors.  The  construction  of 
the  diagram  described  in  this  paper  yields  a  motion-planning  algorithm  for  the 
ladder  which  runs  in  0{n^  log  n  log*/i)  time. 

1.  Introduction 

A  generalized  'Voronoi  diagram'  for  a  line  segment  (ladder')  moviog.  in  2-dimensional  space 
amidst  polygonal  barriers  was  introduced  in  a  companion  paper  by  the  authors  [OSY],  as  a  tool 
for  obtaining  efficient  algorithms  for  planning  obstacle-avoiding  motions  of  the  ladder.  That  paper 
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contains  a  topological  analysis  of  such  diagrams  and  shows  how  to  use  them  for  motion  planning. 
In  the  present  work  we  consider  the  problem  of  how  to  construct  these  diagrams  efficiently.  We 
show  that  if  the  space  in  which  the  ladder  is  confined  to  move  is  bounded  by  obstacles  whose 
boundaries  consist  of  n  straight  segments  and  comers,  then  the  size  of  (a  skeleton  representation 
of)  the  associated  Voronoi  diagram  is  at  most  0(n'  log*n),  and  that  this  diagram  can  be 
constructed  in  time  0(n'  log  n  log'n),   yielding  motion  planning  algorithms  of  similar  time 

complexity.  (As  usual,  log*/i  is  the  height  of  the  smallest  tower  of  exponentials  2-  whose  value 
exceeds  n.) 

The  terminology  used  in  this  paper  was  introduced  in  [OSY],  but  for  convenience  we  review 
it  here.  Let  O  denote  the  free  open  2-dimensional  space  in  which  the  ladder  B  =  PQ  is  confined 
to  move,  let  d  denote  the  length  of  B,  and  let  S  denote  the  set  consisting  of  all  comers  and  (open) 
segments  constituting  the  boundary  of  Cl.  Let  FP  denote  the  3-dimensional  configuration  space  of 
free  placements  of  B  (i.e.  placements  at  which  the  ladder  is  wholly  contained  in  Cl).  Each 
placement  Z  (.  FP  is  specified  by  a  triple  (x^,9),  where  (jc^)  e  ft  is  the  position  of  one  of  the 
endpoints  P  of  B  and  9  is  the  orientation  of  fl;  Z  has  the  property  that  when  B  is  given  this 
placement,  it  lies  wholly  in  H.  By  B(Z)  we  will  denote  the  subset  of  fl  occupied  by  B  when  it  is 
given  the  placement  2.  Note  that  FP  is  an  open  3-dimensional  manifold.  The  space  SFP  of  semi 
free  placements  of  the  ladder  is  defmed  to  consist  of  placements  Z  at  which  B{Z)  is  contained  in 
the  closure  of  Cl.  SFP  is  a  3-dimensional  manifold  with  boundary.  We  reserve  the  terra  "serai- 
free"  to  refer  to  placements  in  SFP-FP. 

The  generalized  Voronoi  diagram  Vor^{S)  for  the  ladder  B  moving  amidst  the  obstacles  of  S 
is  defmed  as  foDows.  For  each  Z  i  FP  let  CUarance{Z)  =  rain  {d(u,v)  :  u  €  B(Z),  v  €  Cf),  and 
let  Near{Z)  be  the  set  of  all  points  in  fy  at  which  the  distance  Clearance  (Z)  is  attained.  Then 
Vor^^S)  is  defined  to  be  the  set  of  ail  placeraents  Z  i  FP  for  which  Near(Z)  is  disconnected. 

It  is  shown  in  [OSY]  that  Vorj(S)  is  a  two-dimensional  subset  of  FP.  We  partiticm  it  bto 
Voronoi  sheets  defmed  as  follows.  For  each  pair  of  separated  obstacles  J;,  J;  €  S  we  defmc  the 
midsurface  M(s~^  ,^2)  to  consist  of  those  free  placements  Z  which  are  simultaneously  nearest  to  j. 
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and  to  Ji,  but  to  no  other  object  in  5.  (Z  is  nearest  to  ^  C  5  if  Near(Z)  intersects  s.)  A  sheet  is 
then  defined  to  be  the  closure  in  SFP  of  a  connected  component  of  such  a  midsurface.  If  one 
assumes  that  the  obstacles  in  5  are  in  "general  position"  (for  a  precise  definition  of  which  see 
[OSY])  then  Vor^(S)  can  be  seen  to  be  eqiial  to  the  intersection  of  FP  with  the  union  of  all  such 
sheets. 

The  boundary  of  a  sheet  is  composed  of  edges,  linked  together  by  (Voronoi)  vertices. 
Roughly  speaking,  a  free  (resp.  semi-free)  edge  is  a  connected  1-dimensional  set  of  free  (resp. 
semi-free)  placements  simultaneously  nearest  to  three  (resp.  two)  distinct  objects  in  S.  Similarly  a 
free  (resp.  semi-free)  vertex  is  a  free  (resp.  semi-free)  placement  simultaneously  nearest  to  four 
(resp.  three)  distinct  objects  in  S  (see  [OSY]  for  refinements  and  special  cases  of  this  definition). 

It  is  shown  in  [OSY]  that  the  boundary  of  each  sheet  2  consists  of  at  most  two  connected 
components,  and  these  two  components  can  be  connected  together  by  a  streamline  consisting  of  all 
placements  Z  €  2  having  a  fixed  canonical  orientation,  as  described  in  [OSY]  (where  such  a 
streamline  is  called  a  core  streamline).  We  will  refer  to  the  union  of  the  boundary  of  2  with  this 
streamline  as  the  (semi-free)  skeleton  of  2.  It  is  shown  in  [OSY]  that  the  free  skeleton  of  2, 
defined  as  the  closure  of  the  intersection  of  its  semi-free  skeleton  with  FP,  is  also  connected. 
Hence,  as  shown  in  [OSY],  the  union  N  of  the  free  skeletons  of  all  sheets  in  Vor^^S)  has  the  same 
connectivity  as  Vor^(S).  We  will  call  N  the  skeleton  of  Vor^(S);  note  that  ^  is  a  one-dimensional 
graph  consisting  of  vertices  connected  by  curves. 

The  main  result  of  [OSY]  is 

Theorem:  Let  Z^,  Zj  €  F/».  There  exists  a  continuous  free  motion  of  B  between  Zj  and  Zj  if  and 
only  if  there  exists  a  continuous  motion  of  B  between  ha(Z{)  and  JmiZ^,  which  passes  only 
through  configiirations  of  N,  where  Im  is  a  certain  retraction-like  mapping  of  FP  onto  N. 

Thus  the  object  that  we  need  to  construct  in  order  to  plan  motions  using  this  theorem  is  the 
1-dimensional  network  S  consisting  of  Voronoi  edges  and  core  streamlines.  Once  N  has  been 
constructed,  we  can  always  seek  a  continuous  path  in  AT  in  time  0(\N\)  by  a  simple  graph  search. 
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This  paper  will  show  that  the  size  oi  N  is  0{n'  log'n),  and  that  N  can  be  constructed  in  time 
0(n-  log  n  log'n).  This  will  be  achieved  by  first  computing  all  the  vertices  of  N,  and  then  by 
showing  how  to  determine  when  two  of  these  vertices  are  adjacent  along  edges  oi  N.  All  this  will 
be  seen  to  lead  to  an  0{n^  log  n  log'n)  motion-planning  algorithm  for  a  ladder  in  polygonal  2- 
space,  thus  substantially  improving  the  0{n^)  algorithm  given  in  [SS].  We  will  also  show  that  FP 
can  have  il(n^  connected  components  in  the  worst-case,  which  indicates  that  our  algorithm  is 
dose  to  optimal  if  we  are  required  to  compute  all  connected  components  of  FP  before  motion- 
planning  actxially  begins. 

In  the  following  analysis  we  find  it  convenient  to  classify  the  vertices  of  N  into  several  types 
and  subtypes  and  to  compute  vertices  of  each  kind  separately  by  a  different  procedure.  A  first 
coarse  classification  of  vertices  is  as  follows: 

(1)  Free  Voronoi  vertices  (or  Voronoi  vertices  for  short).  These  are  free  placements  Z  ii  FP  lot 
which  Near(Z)  meets  the  closures  of  four  distinct  objects  in  5,  and  furtermore,  if  j  is  a  wall  edge 
whose  closure  meets  Near{Z),  then  the  distance  of  B{Z)  from  j  is  the  same  as  the  distance  from 
B(Z)  to  the  infmite  line  containing  s. 

(2)  Semi— free -Voronoi  vertices.  These  are  placements  Z  €  N—FP  for  which  5(Z)  touches  30  at  a 
set  of  points  which  meets  the  closures  of  three  distinct  objects  in  5,  with  the  additional 
requirement  as  in  (1)  for  the  wall  edges  among  these  objects. 

(3)  Pseudo  Voronoi  vertices.  These  are  placements  which  belong  to  the  intersection  of  a  core 
streamline  (see  [OSY])  and  a  Voronoi  edge  in  the  ladder's  configuration  space. 

The  free  Voronoi  vertices  are  hardest  to  compute  efficiently,  so  most  of  this  paper  will  be 
devoted  to  their  calculation.  (In  this  and  subseqtience  development  we  will  only  deal  with  vertices 
dcfmed  by  four  separated  objects.  At  the  end  of  Section  4  we  will  remark  on  how  the  methods 
developed  in  the  paper  extend  to  all  other  cases.)  The  free  vertices  can  be  classified  as  follows  (cf. 
Fig.  1.1). 

Vertices  of  Type  I:  These  are  free  placements  Z  for  which  fl(Z)  is  nearest  to  four  separated  objects 
1],  J;,  jj  and  j^,  such  that  the  distances  from  fl(Z)  to  s.^  and  s^  are  both  attained  at  the  same 
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endpoint  of  B  (two  objects  s^,SivaS  art  said  to  be  sqjarated  from  one  another  if  they  are  neither 
a  segment  and  one  of  its  endpoints  nor  two  endpcints  of  the  same  segment). 

Note  that  in  this  case  we  do  not  further  restrict  the  points  on  fl(Z)  at  which  its  distances  from  jj 
and  s^  are  attained,  so  that  in  particular  these  can  also  coincide  with  the  endpoint  of  B 
simultaneously  nearest  to  s-^  and  t^,  can  lie  along  the  sides  of  B{Z),  etc 

All  Voronoi  vertices  Z  not  falling  into  this  rather  general  case  will  be  referred  to  as  vertices 
of  Type  n.  In  such  configurations  each  endpoint  of  5(7)  can  be  nearest  to  at  most  one  object  in  5, 
and  thus  at  least  two  other  objects  in  S  lying  nearest  to  5(Z)  must  be  nearest  to  interior  points  of 
S(Z).  We  can  rule  out  the  case  of  four  objects  in  5  all  being  nearest  to  interior  points  of  5(Z), 
and  also  the  case  of  three  objects  in  5  all  being  nearest  to  interior  points  of  5(Z)  and  lying  on  the 
same  side  of  5(Z),  since  it  is  easily  seen  that  this  would  violate  our  assumption  that  the  obstacles 
are  in  "general  position".  It  follows  that  all  TVpe  11  vertices  must  fall  into  one  of  the  following 
subcases  (sec  Fig.  1.1):  — 

Vertices  of  Type  H.a:  These  are  free  placements  Z  such  that  5(Z)  is  nearest  to  four  separated 
objects  s^,  S2,  s^  and  s^,  such  that  the  distances  from  B(Z)  to  j^  and  ^2  are  both  attained  at  interior 
points  of  B(Z),  with  both  Sy  and  ^2  ^yuig  on  the  same  side  of  B(Z),  and  such  that  the  distances 
from  B(Z)  to  s^  and  J4  are  attained  at  different  endpoints  of  B. 

Vertices  of  Type  II. b:  These  are  free  placements  Z  such  that  B(Z)  is  nearest  to  four  sqjarated 
objects  Sy,  S2,  53  and  .^4,  such  that  the  distances  from  B(Z)  to  j^,  ^2  and  s^  are  all  attained  at 
interior  points  of  B(Z),  with  s^  and  ^2  lying  on  one  side  of  B(Z),  and  with  J3  lying  on  the  other 
side. 

Vertices  of  Type  II. c  These  are  free  placements  Z  such  that  5(Z)  is  nearest  to  four  separated 
objects  j^,  S2,  J3  and  s^,  such  that  the  distances  from  5(Z)  to  j^  and  J2  ^^  both  attained  at  interior 
points  of  5(Z),  with  j^  lying  on  one  side  of  B  and  52  lying  on  the  other  side,  and  such  that  the 
distances  from  B(Z)  to  jj  and  s^  are  attained  at  different  endpoints  of  B. 

Remark:  In  the  above  classification  we  have  assumed  that  the  four  objects  s^,  .  .  .  ,s^  involved  are 


separated.  This  is  done  to  simplify  the  presentation  of  the  various  algorithms  for  calculating  these 
vertices  in  the  following  sections.  However,  we  will  note  later  that  much  the  same  methods  also 
apply  to  the  cases  where  not  all  these  objects  are  separated. 


Type  n . b 


Fig.  1.1.  Free  Voronoi  vertices 
The  paper  is  organized  as  follows:  Section  2  describes  a  procedure  for  computation  of 
Voronoi  vertices  of  type  I,  Section  3  discusses  the  computation  of  vcrtiocs  of  type  n.a  and  n.b, 
and  Section  4  shows  how  to  find  vertices  of  type  II.c.  A  technique  for  fmding  all  semi-free  vertices 
and  all  pseudo-vertices  is  given  in  Section  5,  and  finally  a  technique  for  constructing  the  edges  of 
the  skeleton  N  is  described  in  Section  6. 

Note  finally  that,  although  our  interest  in  introducing  and  analyzing  generalized  Voronoi 
diagrams  of  this  sort  is  motivated  by  their  application  to  motion-planning,  we  believe  that  the 
concept  is  of  independent  interest. 
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2.  Voronol  Vertices  of  Type  I 

In  this  section  we  only  consider  Voronoi  vertices  of  Type  I,  prove  that  there  are  at  most 
0(n^  log'n)  such  vertices,  and  present  an  algorithm  for  computing  all  these  vertices  in  time 
0(n^  log  n  log'n),  where  as  before  log'n  is  the  smallest  integer  i  such  that  EXPi  >  n,  where 

£X/»i  =  2  and  EXPj+i  =  2  '.  Subsequent  sections  will  consider  other  Jypes  of  free,  semi-free, 
and  pseudo-vertices  and  will  establish  similar  bounds  on  the  numt)er  of  sucfa  vertices  and  on  the 
time  needed  to  calculate  them. 

Voronoi  vertices  of  type  I  are  placements  Z  of  the  ladder  B  for  which  B(Z)  is  nearest  to  four 
objects  Sj,  7=1,  ...  ,4  of  5  and  such  that  there  exist  at  least  two  of  these  objects,  which  for 
definiteness  we  will  call  s^  and  J2»  such  that  d(si,B(Z))  and  d(s2,B(Z))  are  both  attained  at  the 
same  endpoint  P  of  B.  To  each  such  Voronoi  vertex  Z  there  corresponds  an  edge  e=e,  of  the 
standard  Voronoi  diagram  Vofq^S)  of  S  (see  Kirkpatrick  [Ki]  for  a  discussion  of  the  properties  of 
this  standard  diagram  for  sets  consisting  of  line  segments  and  points  in  the  plane).  This  edge 
consists  entirely  of  points  in  the  plane  simultaneously  nearest  to  s^  and  S2,  and  P  must  dearly  lie 
on  tfj  ^ .  See  Fig.  2.1  for  a  description  of  the  possible  forms  of  such  a  standard  Voronoi  edge. 
We  will  refer  to  Z  as  a  vertex  (of  Vor^(S))  generated  by  s^  and  53  and  lying  along  e  =  e,  ^.  [Ki] 
shows  that  the  total  number  of  edges  in  Voro(J)  is  0(ji),  and  our  first  aim  will  be  to  show  that  for 
each  edge  «=«,  ^  in  Votq^S)  there  can  exist  at  most  0{n  log*/i)  vertices  of  VorJ^S)  generated  by 
s-^  and  ^2  and  lying  along  e. 

Thus  let  jj,  ^2  be  given,  and  let  «  be  an  edge  of  Voro(5)  all  of  whose  points  are  equidistant 
from  jj  and  S2-  For  each  X  i  e  define  d{X)  =  d(X,s{)  =  d{X,s^.  Since  jj  and  jj  can  be  either 
edges  or  comers,  we  need  to  consider  three  possible  types  of  edges  e,  for  each  of  which  we  will 
assume  that  j^  and  ^2  are  placed  in  an  appropriate  canonical  position  (see  Fig.  2.1).  The  three 
possible  cases  are  as  follows. 

(CC)  Both  jj  and  ^2  are  wall  comers,  in  wfaidi  case  we  place  s^  and  ^2  on  the  >  axis,  and  «  on  the 
X-axis.    In  this  case  the  origin  M  divides  *  into  two  portions,  along  each  of  which  the  function 
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d(X)  is  monotone. 

(WQ  jj  is  a  wall  edge  and  jj  u  a  wall  comer.  In  this  case  we  place  j^  on  the  j:  axis  and  J2  on  the 
positive  y  axis.  Again  the  point  M  on  the  y  axis  halfway  between  J2  andThe  x-axis  divides  e  into 
two  portions,  along  each  of  which  d{X)  is  monotone. 

(WW)  Both  jj  and  J2  ore  wall  edges  (which  we  assume  to  be  nonparallel).  Here  we  can  assume 
that  the  lines  containing  j^  and  53  meet  at  the  origin,  that  both  s^  and  jj  lie  in  the  right  half-plane 
and  that  the  positive  x-axis  is  the  bisector  of  the  angles  between  them.  In  this  case  d(X)  is 
monotone  along  the  entire  edge  e. 


X 


Fig.  2.1.  Canonical  forms  of  standard  Voronoi  edges. 
The  following  definitions  (cf .  [OSY])  will  be  useful: 

DeflnltioD:  Let  Z  =  (xj,9)  be  a  free  placement  of  the  ladder,  and  let  r  >  0.  The  r-racetrack 
about  Z,  denoted  as  R{Z,r),  is  the  set  of  all  points  whose  distance  from  B{Z)  is  Sr.  The  boundary 
of  this  racetrack  R  consists  of  two  straight  segments  (called  the  sides  oi  R)  parallel  to  B(Z)  at 


distance  r  apart,  and  of  two  sani-drcles  of  radius  r  (called  the  bends  of  R)  about  P  and  Q 
connecting  the  two  sides.  R  is  called  cUar  if  no  object  of  S  meets  the  interior  of  R.  By  R(Z)  we 
denote  the  maximal  clear  racetrack  about  Z. 

Definition:  Let  j  be  an  obstacle  (i.e.  a  wall  or  comer),  and  let  6  be  an  orientation.  By  5,  we 
denote  the  set  {(r  cos9,  r  sin9)  :0:Sr^d},  and  we  define  A(j,9)  to  be  the  Minkowski  difference 

s  -  Bf  =  {x  -  y:x  ^  s,y  (^  B^ 
Let  tQ  be  one  of  the  portions  of  e  along  which  d(X)  is  monotone,  and  consider  the  problem 
of  finding  the  type  I  Voronoi  vertices  (X,Q)  generated  by  s^  and  ^2  for  which  X  lies  on  cq. 

First  note  that  if  X  is  simultaneously  nearest  to  s^,  52,  and  a  third  object  53,  then  X  is  a 
standard  Voronoi  vertex  bounding  e^.  For  each  such  vertex  X  we  can  easily  find  all  orientations  9 
for  which  B(X,Q)  is  simultaneously  nearest  to  the  three  objects  Sj,j=l,  ...  ,3  at  Jf,  and  to  a 
fourth  object  J4  €  5  at  some  other  point  on  it,  as  follows.  Let  r  denote  the  common  value  of  the 
distances  d{X,sj),  7=1,2,3.  Apply  a  standard  technique  whidi  sweeps  the  racetrack  R{{Xfi),r) 
about  X  by  rotating  9  a  full  drcle,  and  also  mnintfliTw  a  list  of  all  objects  in  S  that  intersect  the 
interior  of  the  racetrack.  It  is  dear  that  this  list  has  to  be  updated  at  most  0{n)  times,  and  that 
when  performing  each  of  these  updates  we  can  check  whether  the  racetrack  is  dear,  in  which  case 
a  Voronoi  vertex  of  the  sort  we  seek  will  have  been  foimd.  This  procedure  can  dearly  be 
implemented  to  run  in  time  0(ji  log  n),  and  since  the  total  number  of  vertioes  of  Voro(S)  is  0{n) 
(cf.  [Ki]),  it  follows  that  all  Voronoi  vertices  in  VorJ^S)  of  the  type  just  discussed  can  be  found  in 
time  0{n^  log  n).) 

We  can  therefore  assume  that  X  is  an  interior  point  of  the  standard  Voronoi  edge  ^g.  We  can 
also  assume,  without  loss  of  generality,  that  «o  lies  in  the  right  half-plane.  Let  A/  denote  the 
leftmost  point  on  e^,  i.e.  the  point  X  ^  «q  at  which  d{X)  a  minimal. 

In  all  three  of  the  cases  listed  above,  if  Z  =  (X,9)  is  a  placement  of  the  ladder  B  such  that 
X  i  Cq,  9  e  (-1T/2,  IT/2),  and  d{i-^fi{Z))  =  d{s2fi{Z))  then  9  must  lie  in  an  interval  [a,P], 

defined  by  X,  where  --^  <  arz^  <  ^.  This  interval  is  constant  m  case  (WW),  and  otherwise 
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grcws  as  X  moves  away  from  A/  along  e^. 

Deflnldon:  Let  j  €  5  be  a  third  object.  For  each  6  €  {-v/1,  v/2)  we  denote  by  X/S)  the  position 
X  €  «o  (if  any)  for  which 

d(B(X,9),s)  =  diB(X,e),Sj)  ,  j=l,2.  (1) 

and  for  which  the  distances  from  B(X,9)  to  5;  and  ^2  ^c  attained  at  the  endpoint  P  of  B,  whereas 

the  distance  from  5(^,9)  to  j  is  attained  at  another  point  of  B  (sec  Fig.  2.2). 


Fig.  2.2.  The  critical  position  X,(6). 
Note  for  use  below  that  the  point  X/9)  is  equidistant  from  the  three  sets  s,,  j;,  and  A(j,9). 
Note  also  that  X,(9)  having  the  properties  appearing  in  the  preceding  definition  need  not  always 
exist.  For  example,  the  distance  from  every  point  X  on  «o  to  A(j,9)  may  be  less  than  the  distaiKc 
from  any  such  X  to  jj  and  J2.  or  for  all  such  X  the  distance  from  X  to  A(5,9)  may  be  larger  than 
that  to  J;  and  j;-  In  the  first  case  we  put  X,(9)  =  -»  and  in  the  second  case  we  put  X,(9)  =  +=». 
It  may  also  be  the  case  that  the  domain  of  definition  of  X,(9)  is  not  connected.  However  it  is  plain 
that  this  domain  of  definition  must  consist  of  some  finite  number  of  intervals,  the  total  number  r 
of  which  is  bounded  by  some  constant  independent  of  n  (f  depends  only  on  the  degrees  of  the 
algebraic  equalities  and  ineqxialities  defining  X,(Q)).  Thus  we  can  always  represent  each  X,(9)  as 
the  minimum  of  a  constant  number  of  functions,  eadi  defined  on  some  subinterval  /  of 
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(-11/2,  TtH),  possibly  equal  to  -»  on  /,  and  pxjssibly  equal  to  +"  outside  /. 

Lemma  2.1:  (a)  For  each  s  (.  S  and  each  fixed  9  C  (-11/2,  it/2)  there  exists  at  most  one  point 
Xj(6)  in  gQ  having  the  properties  appearing  in  the  preceding  definition. 

(b)  Let  J  6  5  and  9  €  (-it/2,  it/2)  be  given  and  suppose  that  X^(9)  is  defined.  Then  for  points 
X  i  eQ  left  of  X,(Q)  we  have  d(X,  s-^  <  d{X,L{sfi))  and  conversely  for  points  X  (.  e^  right  of 
XX9)  we  have  d{X,  jj  >  dQCAiifi))- 

Proof:  Let  e  denote  the  locus  of  points  equidistant  from  s^  and  ^2  (so  that  ^g  is  a  subset  of  e). 
Consider  the  function 

<J>(;ir,9)  =  d(X,  jj)  -  d{X,  A(j,9))  ,    X  ^e 
Note  that  ^  is  continuous  on  e.    Let  X  ^  e^ha  a  point  satisfying  (1),  so  that  4>(X,9)  =  0.   We 

shall  show  that  f or  9  C  (-ir/2,  it/2)  the  directional  derivative  of  ^(Xfi)  must  be  positive  at  every 

zero  of  4>(^i9)  as  we  move  right  along  the  curve  e^,  from  this  it  follows  that  <)>(X,9)  has  at  most 

one  zero  X^ifi)  on  e^,  and  that  X,(9)  has  the  property  (b).  To  prove  this  assertion  concerning  the 

derivative  of  ({>,  note  first  of  all  that  the  minimum  distances  of  X  from  j^,  s^  ^^"^  ^i  ~  ^i^fi)  ^^ 

attained  at  points  PiJJ2J>i  of  *i»  •»2»  *3  respectively,  such  that  for  each  j,  the  set  Sj  lies  on  that  side 

of  the  line  perpendicular  to  Xpj  at  pj  which  lies  opposite  to  X.  Hence  it  follows  that  if  X  is  moved 

to  y  by  an  infinitesimal  displacement,  the  distances  of  Y  from  Si,S2t'i  respectively  are  equal  up  to 

first  infinitesimal  order  to  Y's  respective  distances  from  PiJ>2J'i-  From  this  it  follows  that  in  every 

one  of  the  three  cases  shown  in  Fig.  2.1,  the  tangent  line  to  the  curve  cq  at  the  point  X  i  cq 

bisects  the  angle  PiXp2,  and  that  when  given  the  placement  (X,Q)  the  ladder  (which  it  is 

convenient  to  regard  as  the  directed  segment  PQ  from  P  to  Q)  points  out  of  this  angle. 

It  is  easiest  to  follow  the  remainder  of  our  argument  by  thinking  of  the  situation  in 
coordinates  like  those  used  in  the  left  hand  portion  of  Fig.  2.1,  in  which  p^  and  pj  lie  on  the  y-axis 
symmetrical  to  one  another,  X  lies  on  the  positive  x-axis,  and  the  tangent  to  «q  at  X  is  the  x-axis. 
By  the  preceding  remarks,  this  involves  no  loss  of  generality.  (These  are  not  quite  the  same  as  the 
coordinates  used  in  the  middle  and  right  portions  of  Fig.  2.1.) 
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Let  r  be  the  cxjinmon  value  of  the  distances  d(X,s{),  d{X,s-^,  d(X,A(s,Q))  at  the  zero  X  €  *o  of 
<j>(Ar ,6).  Then  pj  lies  on  the  aide  C  of  radius  r  about  X.  We  shall  show  that  in  fact  p^  must  lie  on 
the  larger  arc  A  of  the  two  arcs  a,  A  into  which  the  points  pj,  P2  divide  this  circle.  To  sec  this, 
observe  that  because  the  cndpoint  /»  is  the  point  of  the  ladder  closest  to  both  pj  and  to  p^,  0  is 
confined  to  the  angle  formed  by  the  rising  tangent  vector  to  C  at  p^  and  the  falling  tangent  vector 
to  C  at  />2-  Furthermore,  p^  cannot  belong  to  j,  for  that  would  make  X  eqddistant  from  the  three 
objects  J,  Jj,  and  Ji.  so  that  X  would  be  a  vertex  of  Vor^^S),  contrary  to  assumption.  Hence  there 
must  exist  a  point  W  ^  s  such  that  W  =^  pj  and  such  that  the  segment  pyW  has  orientation  6.  But 
since  this  segment  is  wholly  contained  in  A(j,9),  it  cannot  pass  through  the  interior  of  C.  If  the 
point  pj  were  on  the  left  hand  arc  a  of  the  drde,  any  oriented  segment  p^W  pointing  right  at  an 
angle  lying  between  the  tangent  vectors  to  a  at  its  endpoints  would  dearly  have  to  pass  through 
the  interior  of  C.  Hence  p^  must  lie  on  A,  as  asserted.  (Note  however  that  it  is  possible  for  p^  to 
coindde  with  one  of  p^  or  pj,  say  for  definiteness  with  p^,  but  then  6  must  be  in  the  direction  of 
the  tangent  to  C  at  p^.  It  is  also  easy  to  see  that  both  s^  and  s  must  lie  on  that  side  of  the  line 
through  pi  at  orientation  9  which  lies  opposite  to  X.  But  then,  as.suming  that  the  objects  in  5  are 
in  "general  position",  it  is  easily  seen  that  neither  jj  nor  s  can  be  a  wall  edge,  for  otherwise  we 
would  have  collinearity  between  a  wall  edge  and  a  non-adjacent  wall  comer,  contrary  to 
assumption.  Hence  in  such  a  coincidence  both  s^  and  s  must  be  corners  and  the  directed  segment 
s^s  must  have  orientation  6.) 

We  therefore  condude  that  the  scalar  product  of  the  vector  Xp^  with  the  unit  tangent  t  to  e^ 
at  X  (which  in  each  of  the  cases  shown  in  Fig.  2. 1  we  take  to  point  to  the  right)  is  strictly  larger 
than  the  scalar  product  of  t  with  either  Xp^  or  Xp^,  these  two  latter  products  being  equal  (this  is 
true  even  in  the  case  where  P3  coinddcs  with  either  p^  or  P2,  as  can  be  verified  directly).  Oir 
original  assertion  about  the  directional  derivative  of  4)(X,6)  along  tQ,  and  with  this  the  present 
lemma,  follow  immediately.  Q.E.D. 

Each  placement  (X,(9),9)  in  which  the  ladder  lies  doscr  to  jj,  jj,  and  s  than  to  any  other 
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obstacle  lies  along  a  Voronoi  edge  in  VorJ^S).  For  fixed  6,  Lemma  2.1  implies  that  for  points 
X  i  eQ  whose  first  coordinate  x  is  greater  than  the  first  coordinate  0:^(8)  of  X,(9),  the  ladder  will 
be  closer  to  s  than  to  s-^  and  J2-  "^lus  in  seeking  for  Voronoi  vertices  generated  by  jj  and  jj  we 
need  only  be  concerned  with  placements  (Xj(9),8)  for  which  Xj(9)  =  minz,.(9),  where  s'  ranges 

over  the  set  S'  of  all  remaining  obstacles.  Moreover,  such  a  placement  is  a  Voronoi  vertex  in 
Vor^(S)  precisely  when  the  ladder  lies  simultaneously  nearest  to  four  obstacles,  i.e.  when  there 
exist  two  objects  s,  s*  such  that 

'  1  i  S 

Thus  to  find  the  desired  vertices  of  Vor^(S)  we  only  have  to  compute  the  function 

A/(e)  =    mn  x,.(9)  , 
t  (  s 

and  to  find  all  those  9's  at  which  this  minimum  is  simultaneously  attained  by  two  functions  x^(9), 
x.(9). 

An  upper  bound  on  the  number  of  such  9's  can  be  derived  by  using  an  interesting 
combinatorial  result  for  sequences  that  do  not  contain  certain  forbidden  subsequences,  as  derived 
by  Szemeredi  [Sz],  Davenport  and  Schinzel  pS],[Da],  and  recently  by  Atallah  [At]  and  by  Hart 
and  Sharir  [HS].  More  specifically,  let  s  ^  S'  and  let  /  be  the  set  of  9  €  (-ir/2,  it/2)  at  which 
^si^)  '^  ■^°°  (i-s-  is  either  defined  or  is  -»).  As  noted  we  can  assume  without  loss  of  generality 
that  /  is  a  connected  subinterval  of  (-17/2,  tt/I).  Let  9  be  the  left  (resp.  right)  endpoint  of  I. 
Extend  the  graph  of  Xj(jd)  to  the  left  (resp.  right)  of  9  as  a  decreasing  (resp.  increasing)  straight 
line  with  a  very  steep  slope,  thus  making  it  continuous  over  the  entire  interval  (-11/2,  'it/2). 
Similarly,  if  /  is  a  subinterval  of  /  at  which  Xj(9)  =  -",  modify  the  graph  of  Xj(9)  on  /  by 

rqjladng  it  by    a  very  steeply  decreasing  line  segment  followed  by  a  very  steeply  increasing 

ft 
segment.   It  is  plain  from  the  fact  that  each  of  the  functions  x,  is  piecewise  algebraic  (in  tan—) 

and  of  a  degree  known  a  priori  that  any  pair  of  the  functions  Xj(6)  as  thus  extended  intersect  in  at 
most  k  points,  where  ifc  is  a  constant  depending  only  on  the  degree  of  the  functions  x^,  but  not  on 
n.  It  is  then  easy  to  see  (cf.  [DS],  [At])  that  the  number  K  of  values  9  at  which  the  minimum 
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A/(9)  is  simultaneously  attained  by  two  functions  x/8)  and  x^.(9)  is  at  most  the  length  \{n,k)  of 
the  longest  sequence  a  each  of  whose  components  belongs  to  {1,  .  .  .  ./>},  which  satisfies  the 
following  two  conditions: 

(i)  For  each  i,  cr^  =^  a,+y 

(ii)    There    do    not    exist    indices    «i  <  «2  <   •  •  "   <  «i+2    **"*    *^^    ^i.  -  '^ij-  ^i,~    '  '  '  • 

a,  =  a,  =  a,  =    •  •  •  i  and  a,  ^  a,. 

To  see  this  partition  (-tt/2,  -n/l)  into  a  minimal  number  of  subintervals  /j I/^  such  that  on 

each  such  interval  /,  the  minimum  A/(9)  is  attained  just  by  a  single  function  x^,,  and  define  ct  to  be 

(jf  ^'....  ,  Jjf).  It  is  then  easy  to  see  that  a  has  properties  (i)  and  (ii),  since  between  any  pair  of 
points  at  one  of  which  x,  <  x .  and  at  the  other  of  which  x,  >  x .,  there  must  exist  a  point  at 
which  these  two  functions  intersect. 

The  combinatorial  results  of  Szemeredi,  Davenport  and  Schinzel,  Hart  and  Sharir,  and 
AtaUah  state  that  \{n,l)  =  In  -  I  pS],[At],  X(n,3)  =  0{na{n))  [HS],  and  \{n,k)  =  0{n  log'n) 
for  each  Jl:  >  3  (with  the  constant  of  proportionality  depending  on  k)  [Sz].  (Here  a(n)  is  the 
functional  inverse  of  Ackermann's  function  which  grows  extremely  slowly;  see  [HS]  for  more 
details).  It  follows  that  K  =  0{n  log'n).  Thus,  since  there  can  be  at  most  0(n)  edges  ^q  in  the 
standard  Voronoi  diagram  VorQliS),  it  follows  that 

Theorem  2.2:  There  are  0(n'  log'n)  Voronoi  vertices  of  type  I. 

Having  now  bounded  the  number  of  Type  I  Voronoi  vertices  which  can  exist,  we  go  on  to 
show  that  all  these  vertices  can  be  found  in  time  0{n-  log  n  log'n)  using  the  following 
straightforward  technique.  We  first  compute  Vor^iS)  in  time  0{n  logn),  using  Kirkpatrick's 
algorithm  [Ki],  or  in  time  0{n  log-n),  using  the  algorithm  in  [LD].  For  each  of  the  0{n)  edges 
*o  =  *o(ji,j2)  of  this  diagram  we  then  compute  the  lower  envelope  A/(9)  of  the  functions  x/9) 
defined  above,  using  the  following  divide  and  conquer  approach. 
(I)  Partition  the  objects  in  S  (except  j^,  jj)  into  two  equal  sized  sets  S^  and  Si- 
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(n)  Recursively  compute  the  lower  envelope  A/(  of  the  functions  x/9)  with  j  C  5/,  for  «  =  1,2. 
Each  such  computation  implicitly  partitions  (-it/2,  it/2)  into  subintervals  such  that  A/^  coincides 
with  a  fixed  function  x^(6)  throughout  each  such  subinterval,  and  assigns  a  particular  minimizing 
function  Xj(9)  to  each  such  subinterval. 

(in)  Compute  A/(e)  =  min  (A/i(9),  A/2(9)).  This  is  done  by  merging  the  two  partitions  of 
(-it/2,  ir/2)  corresponding  to  A/j  and  M2,  and  then  by  scanning  (-ir/2,  it/2)  from  left  to  right, 
computing  the  minimum  M  within  each  of  the  refined  subintervals,  each  such  subccnnputation 
requiring  only  constant  time. 

To  estimate  the  time  cost  of  this  procedure,  let  T(n)  denote  the  complexity  of  the  above 
algorithm  in  the  case  when  the  number  of  objects  s  (i.e.  functions  ^^(9))  being  processed  is  n. 
Then  it  is  plain  that 

Tin)  =  2T(n/2)  +  OiK(n)) 
M^ere  K{n)  is  the  maximal  number  of  changes  in  the  corresponding  lower  envelope  A/.  As  shown 

above  we  have  K(n)  =  0(n  log'/i),  from  which  it  follows  easily  that 

T(n)  =  0(n  log  n  log'n)   . 
This  proves 

Theorem  2.3:  The  collection  of  all  Type  I  Voronoi  vertices  can  be  found  in  total  time 
0(n^  log  n  log'n). 

Remark:  The  bounds  established  in  Theorems  2.2  and  2.3  depend  on  the  maximal  number  /  of 
mtersections  of  any  pair  of  the  functions  X,(9).  If  we  knew  that  /  =  2  the  number  K(n)  of 
changes  in  the  lower  envelope  of  the  jr,(9)'s  could  be  at  most  0(n),  which  would  imply  that  the 
total  number  of  type  I  vertices  must  be  0(n^  and  that  all  these  vertices  can  found  in  time 
0(n^  log  n).  Similarly,  if  we  knew  that  /  =  3  the  total  number  of  type  I  Voronoi  vertices  would 
be  at  most  0(n^a{n)),  and  all  these  vertices  could  be  foimd  in  time  0(n^a(n)  log  n).  On  the 
other  hand,  even  if  /  a  4,  the  slightly  weaker  bounds  given  above  still  apply.  The  actual 
determination  of  /  is  non-combinatorial,  and  involves  only  elementary  analytic  geometry.  We 
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leave  precise  determination  of  /  as  an  open  problem;  however,  initial  examination  of  this  problem 
leads  us  to  conjecture  that  7^3. 

3.  Computation  of  Voronoi  Vertices  of  Types  Il.a  and  n.b 

Once  we  have  found  all  the  type  I  Voronoi  vertices  of  Vor^(S)  we  can  use  this  information 
to  compute  vertices  of  all  the  other  types  efficiently. 
We  consider  vertices  of  types  n.a  and  n.b  first. 

Theorem  3.1:  There  are  at  most  0(n^  log'n)  Voronoi  vertices  of  type  n.a  and  n.b,  all  of  which 
can  be  computed  in  time  0(n^og  n  \og'n). 

Proof:  Let  Jj  and  ^2  be  two  comers  in  5.  The  Voronoi  vertices  Z  of  type  n.a  or  H.b  that  are 
defined  by  s^  and  ^2  are  those  placements  of  the  ladder  at  which  B(Z)  is  nearest  to  jj,  Ji  aod  two 
other  objects  of  S,  such  that  the  distances  d{s^,B(Z))  and  d(s2,B(Z))  are  attained  at  interior  points 
of  B(Z),  and  both  s^  and  ^2  lie  on  the  same  side  of  B  (these  placement  are  called  broadside 
placements  in  the  first  part  of  this  paper  [OSY]).  Call  such  a  Z  a  type  n  vertex  generated  by  jj 
and  J2  (again  we  note  that  since  Jj  and  J2  are  assumed  to  be  separated  obstacles  they  both  must  be 
comers  rather  than  line  segments).  For  each  such  Z,  the  line  segment  fl(Z)  must  be  parallel  to  the 
line  /  connecting  j;  and  ^2-  Without  loss  of  generality,  we  can  take  /  to  be  the  x-axis,  with  s.  lying 
left  of  $2,  see  Fig.  3.1.  Fix  the  distaiKC  f>0  of  B  from  /  and  suppose  that  5  lies  on  a  given  side  of 
/  (say  the  upper  half  plane).  Without  loss  of  generality,  assume  also  that  the  directed  segment  PQ 
points  in  the  same  direction  as  the  directed  segment  JjJi-  Consider  the  set  A{t)  of  all  x  sudi  that 
in  placement  Y  =  (x,f,0)  of  B,  R{y,t)  is  touched  by  both  jj  and  s^  and  is  dear.  By  varying  x  the 
racetraclc  can  be  pushed  left  until  the  center  of  its  right-hand  bend  lies  over  s-,,  or  right  until  the 
center  of  its  left-hand  bend  lies  over  Jj.  If  within  this  range  of  possible  positions  a  point  p  of  some 
other  object  encounters  the  left-hand  (rcsp.  the  right-hand)  bend,  it  is  dearly  impossible  to  push 
the  racetrack  so  far  left  (resp.  right)  that  p  will  leave  the  racetrack.  Thus  A{t)  is  either  empty  or  a 
dosed  interval  [/(0.''(0];  furthermore,  since  /?((x,r,0),  t)  increases  steadily  with  f,  /(f)  is 
nondccreasing  and  r{t)  is  nonincreasing  with  t.    Let  tQ  be  the  maximal  t  iot  which  A{t)  is 
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nonempty,  and  let  Xq  =  /(fo),  x^  =  r(t(^.  U  Xq  =  x^  then  either  (xQ.rg.O)  is  a  placement  at  which 
the  racetrack  touches  s^,  s^t  ^^  two  other  objects  (one  on  its  left-hand  bend  and  the  other  on  its 
right-hand  bend),  i.e.  is  a  Voronoi  vertex  of  type  n.a,  or  is  a  placement  at  which  one  end,  say  the 
left  end,  of  the  ladder  lies  above  i-^  and  the  left  bend  of  the  racetrack  touches  a  third  object  of  S. 
Otherwise,  \i  xq^  x^,  then  since  /q  cannot  be  inaeased  without  A(f)  becoming  empty,  it  must  be 
impossible  to  slide  the  racetrack  along  j^  and  jj  to  move  it  out  of  contact  with  the  object  J3  that  it 
touches  along  its  left-hand  (resp.  right-hand)  bend  without  remaining  in  contact  with  some  fourth 
object  s.  Hence  s  must  touch  the  racetrack  on  its  upper  side  for  all  positions  (x,/q,0)  between 
(xfl.ro.O)  and  (oci.fQ.O).  In  this  case  {x^j^fl)  (resp.  (xi,ro,0))  must  be  dther  a  type  n.b  vertex  or  a 
placement  in  uliich  one  end  of  the  ladder  lies  directly  above  5^  or  ^2  while  the  upper  side  of  the 
corresponding  racetrack  touches  s. 

Conversely,  let  Z  =  (x,r,0)  be  a  type  11  vertex  generated  by  s^,  sj.  Then  dther  f  =  fg  or 
t  <  tQ.  If  t  =  (q  then  we  must  have  Z  =  (xQ,tQ,0)  or  Z  =  (xift^fl)  (since  if  Xg  <  x  <  x^  then  the 
upper  side  of  the  racetrack  R(Z)  must  be  touched  by  two  additional  comers  sj,  s^  of  S,  a 
possibility  wiiich  is  ruled  out  by  oiir  assumption  that  the  objects  in  5  are  in  "general  position"  (see 
also  [OSY])).  Next  consider  the  case  in  which  r  <  rg,  so  that  A(t)  has  positive  length.  In  this  case 
the  two  additional  objects  s^,  J4  closest  to  5(Z)  cannot  touch  both  bends  of  R(Z),  nor  can  they 
both  touch  the  same  bend  of  R(Z),  since  this  would  make  Z  is  a  type  I  Voronoi  vertex,  contrary 
to  assumption.  Hence  one  of  these  objects  touches  a  bend  of  R(Z)  and  the  other  touches  its  upper 
side.  Suppose  for  defmiteness  that  jj  touches  the  left-hand  bend  of  R(Z)  and  that  s^  touches  its 
upper  side.  If  s^  lies  to  the  right  of  jj,  then  it  is  dear  that  A^t")  is  empty  for  r*  >  r,  and  thus  s^ 
must  lie  to  the  left  of  s^.  We  thus  conclude  that  every  type  n  vertex  generated  by  jj,  J2  niust 
have  one  of  the  forms  described  by  conditions  (a.2),  (b.O),  and  (b.l)  in  the  following  list  of 
'exceptional'  conditions: 

(a.l)  Two  distinct  objects  J3,  s^  ^  S  touch  the  left  bend  of  the  r-raoetrack  R  about  the 
placement  (l(t),t,0)  of  B.  (In  this  case  the  object  a(r')  touched  in  this  placement  can 
change  as  we  pass  through  t.) 
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(a.2)  One  object  s^  (.  S  touches  the  upper  side  of  R  to  the  left  of  jj  and  another  object  s^  i  S 
touches  the  left  bend  of  R.  (In  this  case  /(r')  will  decrease  discontinuously  as  we  pass 
through  t.  Note  also  that  in  this  case  Z  =  (/(f), f,0)  is  a  Voronoi  vertex  of  type  Il.b.) 

(a.3)  The  right  hand  end  of  the  ladder  lies  directly  above  jj  and  another  object  s  i  S  either 
touches  the  upper  side  of  the  racetrack  R  to  the  left  of  s^,  or  touches  the  left  bend  of  R. 
(These  situations  resemble  cases  (a.2)  and  (a.l)  respectively,  differing  oiJy  in  the  fact 
that  f or  r*  <  f  leftward  motion  of  the  racetrack  is  limited  not  by  the  presence  of  an 
object  which  it  contacts  but  by  the  fact  that  the  right  end  of  the  ladder  lies  over  J2) 

Similar  and  symmetric  conditions  govern  the  oocurrenoe  of  discontinuous  changes  in  r(f)  and  the 
object  3(r)  which  the  racetrack  touches  when  B  is  given  placement  (r{t),t,0). 

The  following  additional  conditions  describe  the  ezcqjtional  cases  whidi  arise  at  r  =  fg,  which  is 

the  maximal  t  for  which  A  (r)  is  not  empty: 

(b.O)       A(tf)  is  a  singleton  and  Z  =  ('('o)='"(f(j),ro,0)  is  a  Voronoi  vertex  of  Type  n.a. 

(b.l)       A(tf^  is  not  a  singleton  and  there  exists  a  comer  j  €  5  siu±  that  for  every  x  in 

['(fo)>''(^o)]i  ^((^.'o>0))  is  cl^^  ^'^  touches  j  on  its  upper  side. 
(b.2)       A(fQ)  is  a  singleton  and  Z  =  (/('q)~''('q)>^0'^)  •*  "  placement  at  which  either  the  left  end 

of  B  lies  directly  above  s^  and  some  third  object  j  C  5  touches  the  left  bend  of  the 

racetrack  /?(Z),  or  a  similar  and  symmetric  condition  holds  for  the  right  end  of  B  and  the 

left  bend  of  ;?(Z). 
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Fig.  3.1.  Critical  heights  in  the  calculation  of  vertices  of  type  n.a  and  Il.b 
To  estimate  the  number  of  these  vertices  and  placements,  and  to  find  them  all  in  an  efficient 
manner,  we  can  proceed  as  follows.  For  each  OsfSfQ,  let  a(r)  (resp.  3(f))  denote  the  set  of 
objects  in  S  which  lie  on  the  left  (resp.  right)  bend  of  the  f-racetrack  about  (/(r),f,0)  (resp. 
(r(f),r,0)).  Note  that  these  sets  may  be  empty  for  small  (possibly  all)  values  of  r.  In  such  cases 
cither  /(f)  coincides  with  jj  -  {dfS)  or  r(f)  coincides  with  s^.  Starting  at  f  =  0,  one  can  gradually 
increase  f  and  track  the  values  of  the  functions  /(f),  rif),  a(f),  P(f).   The  crudal  fact  is  that  as  f 
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varies  the  function  /(f)  is  cxmtinuous  and  a(f)  is  constant  except  for  "critical"  heights  f  at  which 
one  of  the  conditions  appearing  in  the  preceding  list  holds. 

The  observations  just  maHi*.  justify  the  following  procedure.  We  find  all  the  "critical"  heights 
defined  by  conditions  (a.l)-(a.3)  and  process  them  in  increasing  order,  updating  a  and  3  each 
time  f  crosses  a  critical  height,  and  checking  whether  A(i)  can  shrink  to  a  single  point  within  each 
interval  of  noncritical  heights  (since  we  know  a  and  p  for  this  interval,  these  checks  can  be 
accomplished  in  constant  time  as  will  be  shown  in  more  detail  below).  Processing  continues  in  this 
manner  until  either  A(f)  shrinks  to  a  point  (so  that  one  of  conditions  (b.O),  (b.2)  holds),  or  a 
height  fg  satisfying  (b.l)  is  encountered.  Cases  in  which  A(r)  does  shrink  to  a  point  may  be 
Voronoi  vertices  of  type  n.a;  if  if  condition  (b.l)  is  met,  the  two  endpoints  of  A^Iq)  may  be 
Voronoi  vertices  of  type  n.b. 

To  see  that  this  procedure,  which  must  be  applied  to  each  pair  of  comers  in  5,  is  as  efficient 
as  claimed,  we  have  to  bound  the  sum  of  the  total  number  of  critical  heights  over  all  pairs  jj,  ^2  of 
comers  in  S,  and  describe  efficient  techniques  for  finding  these  heights. 

Note  first  that  for  each  t  at  which  condition  (a.l)  holds,  the  placement  (/(f),r,0)  is  a  Voronoi 
vertex  of  type  I.  Thus,  by  the  analysis  presented  in  the  preceding  section,  the  total  number  of  such 
placements,  summed  over  all  pairs  {si,s^,  is  at  most  0(n-  log'n),  and  all  such  placements,  hence 
all  critical  heights  determined  by  condition  (a.1),  can  be  found  in  time  0(n-^og  n  log'n)  using  the 
procedure  already  outlined. 

The  sets  of  critical  heights  satisfying  respectively  conditions  (a.2),  (a.3),  and  (b.l)  can  be 
calculated  by  a  somewhat  different  procedure,  which  in  fact  need  not  yield  precisely  the  set  of 
critical  heights,  but  can  just  as  well  be  allowed  to  produce  a  somewhat  larger  superset,  provided 
that  this  superset  contains  the  set  of  all  critical  heights,  and  provided  that  the  superset  generated  is 
not  too  numerous.  This  is  because  all  spurious  heights  can  simply  be  thrown  away.  (We  emphasize 
again  that  in  the  discussion  which  follows  we  need  concern  ourselves  only  with  configurations 
satisfying  (a.2),  (a.3),  and  (b.l),  but  neither  (a.l)  (since  these  configiu-ations  are  already 
available)  nor  (b.O)  and  (b.2)  (since  these  are  foxmd  "on  the  fly"  by  the  scanning  procedure 
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descxibed  below,  namely  wben  the  interval  A(()  a  seen  to  have  shrunk  to  a  single  point.)  The 
procedures  that  we  are  about  to  present  will  also  be  used  in  calculation  of  Voronoi  vertices  of 
other  types. 

To  find  supersets  of  the  sets  of  placements  Z  satisfying  one  of  the  conditions  (a.2),  (a.3), 
(b.l)  we  consider  the  following  three  geometric  problems. 

Problem  1:  Qven  a  set  5  of  n  points  in  the  plane,  find  all  triples  (A,B,C)  of  points  in  S  having 
the  following  property:  Let  D  be  the  disc  whose  diameter  is  AB.  Then  D  is  required  to  contain  C, 
wWle  the  portion  of  D  bounded  by  the  diameter  AB  and  by  the  line  passing  through  A  and  C  is 
required  to  contain  no  point  of  5  in  its  interior  (see  Fig.  3.2.) 


Fig.  3.2.  A  solution  triple  to  Problem  1 
Problem  2:  Given  a  set  5  of  n  points  in  the  plane,  find  all  triples  (A^,C)  of  pomts  in  S  satisfying 
the  following  conditions: 

(i)  The  dosed  rectangle  ABXY  whose  side  xy  lies  on  the  line  containing  C  and  parallel  to  AS  is 
such  that  XY  contains  C. 

(ii)  The  rectangle  ABXY  and  the  two  adjacent  semidrdes  of  diametera  AY  and  BX  respectively 
contain  no  point  of  5  in  their  interior. 
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(iii)  C  lies  in  the  half-space  to  the  left  of  the  directed  line  containing  the  directed  segment  AS. 
(iv)  The  length  of  AB  is  ^  d. 


Fig.  3.3.  A  solution  triple  to  Problem  2 
Probtem  3:  Given  a  set  S  of  n  points  and  open  line  segments,  and  a  fixed  length  d,  find  all  triple* 
(A,B,C)  where  A  and  B  are  points  in  5  and  C  is  a  point  or  a  segment  in  S  having  the  following 
properties:  Let  AXYZ  be  a  rectangle  svich  that  (i)  the  side  AY  is  of  length  d  and  contains  fl ;  (ii)  the 
midpoint  W  of  XY  is  nearer  to  C  than  to  any  other  object  in  S,  and  the  distance  from  W  to  C  is 
equal  to  WX;  (iii)  the  interior  of  AXYZ  meets  no  object  of  S  (see  Fig.  3.4.) 


Fig.  3.4.  A  solution  triple  to  Problem  3. 


-23- 

The  solution  sets  of  these  three  problems  are  related  in  the  following  way  to  the  set  of 
ladder  placements  satisfying  conditions  (a.2),  (a.3),  or  (b.l)  are  as  follows.  Consider  a  placanent 
satisfying  condition  (b.l),  i.e.  an  empty  racetrack  whose  bottom  side  passes  through  two  points 
s^,  S2  and  whose  top  passes  through  a  third  point  s.  Suppose  for  definiteness  that  the  top  and 
bottom  of  the  racetrack  are  horizontal.  Then  the  projection  of  s  onto  the  line  j^jj  ^^  either 
between  jj  and  ^2.  or  to  the  left  of  j^,  or  to  the  right  of  J2-  First  suppose  that  s  lies  left  of  jj,  and 
draw  the  semicircle  on  the  line  ss2,  as  in  the  following  figure. 


Fig.  3.5.  A  aitical  racetrack  and  the  corresponding  semidrcle. 
It  is  dear  that  both  the  triangle  sxs2  in  Fig.  3.5  and  the  part  of  the  semicircle  sxs2  lying  above  the 

line  51^2  ^^  interior  to  the  racetrack.  Thus  (j,J2,Ji)  constitutes  a  solution  of  Problem  1.  Similarly, 

if  in  condition  (a.2)  we  take  J3  for  j,  it  is  dear  for  the  same  reason  given  above  that  every 

quadruple  satisfying  (a.2)  contains  a  triple  satisfying  Problem  1. 

On  the  other  hand,  if  in  Fig.  3.5  s  lies  vertically  above  the  segment  JiJ2>  then  {s-^^,S2^)  plainly 
constitutes  a  solution  of  Problem  2. 

Finally,  consider  a  placement  satisfying  condition  (a.3).  If  the  point  s  of  condition  (a.3)  touches 
the  racetrack  on  its  upper  side,  then,  as  in  condition  (a.2)  discussed  above,  the  triple  ^21  '>  '1  u  a 
solution  to  Problem  1.  Otherwise,  if  we  let  the  points  jj,  52  of  (a-3)  be  respectively  the  B,  A  of 
Problem  3  and  let  the  point  s  appearing  in  (a.3)  be  the  C  of  Problem  3,  it  is  dear  that  every  such 
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placement  generates  a  solution  of  Problem  3. 

It  is  plain  that  the  number  of  solution  triples  to  each  of  the  three  preceding  problems  is 
0(n^  because  for  each  of  these  problems  and  for  each  pair  A,  B  (^  S  there  exist  at  most  0(1) 
objects  C  satisfying  the  condition  of  the  problem.  On  the  other  hand,  a  straightforward  review  of 
the  simple  constructions  used  just  above  to  relate  placements  satisfying  conditions  (a.2),  (a.3), 
(b.l)  to  solutions  to  Problems  1-3  makes  it  dear  that  for  each  solution  to  one  of  the  problems  1-3 
there  are  at  most  two  placements  satisfying  a  corresponding  condidon  (a.2),  (a.3),  or  (b.l).  This 
implies  that  the  total  number  of  critical  heights  defined  by  conditions  (a.l)-(a.3),  (b.l),  summed 
over  all  pairs  of  comers  s^,  ^2  ^  -5.  is  0(n^  if  we  except  condition  (a.l),  and  without  this 
exception  is  0(n^  log'n). 

We  will  find  all  solutions  to  Problems  2  and  3,  but  will  seethat  to  do  this  is  not  trivial.  In 
fact,  the  methods  to  be  described  below  will  not  yield  exacdy  the  collection  of  all  solution  triples 
to  these  problems,  but  rather  supersets  of  these  collections.  Nevertheless,  each  of  these  supersets 
will  still  be  of  size  0{n^,  and  will  be  calculated  in  time  0(n-^og  n).  The  additional  'spurious' 
heights  present  in  these  supersets  wiH  present  no  problem,  because  when  we  encounter  such  a 
height  we  can  easily  check  whether  it  is  truly  a  critical  height,  and  if  not  simply  discard  it. 

The  following  lemmas  state  that  a  solution  set  to  Problem  1  and  supersets  of  the  solution 
sets  to  Problems  2  and  3  can  be  found  efficiently.  We  will  first  conclude  the  proof  of  Theorem  3.1 
using  these  lemmas,  and  then  gjve  the  proofs  of  these  lemmas. 

Lemma  3.2:  The  collection  of  all  solution  triples  to  Problem  1  can  be  found  in  time  0{n^og  n). 

Lemma  3.3:  A  collection  of  0(n^  triples  containing  all  solutions  of  problem  2  can  be  calculated  in 
time  0(n-^og  n). 

Lemma  3.4:  A  set  of  size  0{n^  which  tnrludfs  all  solution  triples  to  problem  3  can  be  calculated 
in  time  C>(n-^og  n). 

To  condude  the  jjroof  of  Theorem  3.1,  note  first  that  the  union  of  the  solution  sets  to 
Problems  1-3  is  a  superset  of  the  set  of  all  critical  heights  for  a  given  pair  of  points  j,,  5i  €  S.  Let 
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N{s^,s^  be  the  number  of  these  hdghts.  By  preceding  argiiments  2  f^i^i>^2)  -  0(n^og*n). 

Thus  we  can  enumerate  and  sort  all  heights  that  need  to  be  examined  to  find  the  Voronoi  vertices 
of  Types  n.a  and  n.b,  and  this  can  be  done  in  time 

0{n^og  n)  +  2  ^(•»i.J2)  logA^(Ji,J2)  =  OCn^log  n  log'n)  . 
Q.E.D. 

Remark:  It  can  be  shown  that  the  number  of  Voronoi  vertices  of  type  n.a  is  at  most  0{n^ 
(although  it  still  take  0(n^og  n  log'n)  time  to  calculate  than,  as  described  above).  This  is 
because  for  each  pair  s^  and  ^2  of  comers  in  S  there  can  exist  at  most  one  such  vertex  generated 
by  Si  and  s^t  ^i^^^  i^  ^^^^  occur  only  at  the  maximal  height  fg  for  which  A{t^  is  still  nonempty,  and 
that  only  on  one  side  of  the  line  passing  through  s^  and  ^2  (assuming  that  these  comers  are  not 
isolated).  We  do  not  know  whether  the  number  of  type  n.b  Voronoi  vertices  is  also  0(n^, 
because  there  can  exist  several  such  vertices  generated  by  the  same  Si  and  s^,  ocoiring  at  several 
smaller  critical  hdghts. 

Proof  of  Lemma  3.2:  The  following  procedure  will  find  all  solution  triples  {A,B,C)  of  Problem  1 
which  are  such  that  B  lies  to  the  right  of  AC.  A  symmetric  procedure  can  then  be  used  to  find  the 
remaining  triples.  For  each  point  A  i  S,  apply  the  following  procedure:  Sort  all  other  points  in  S 
by  their  angular  orientations  about  A.  Then  iterate  through  the  sorted  list  of  points  in  clockwise 
order.  While  iterating,  maintain  an  auxiliary  ordered  list  of  "active"  points.  Initially,  the  first  point 
processed  is  put  on  the  active  list.  Suppose  that  the  algorithm  is  currently  processing  a  point  X, 
and  let  y  be  the  last  element  on  the  active  list.  If  the  angle  AYX  is  acute,  then  evidently  the  triple 
{A^,Y)  is  not  a  solution  triple  to  Problem  1  (because  Y  must  be  outside  the  circle  whose  diameter 
is  AX).  Moreover,  it  is  also  easy  to  see  that  no  point  Z  lying  further  in  the  clockwise  direction 
from  X  can  define  a  triple  {A^,Y)  which  is  a  solution  to  Problem  1  (because  for  all  such  Z  the 
point  X  must  lie  in  the  portion  of  the  disc  with  diameter  AZ  which  is  bounded  by  AZ  and  by  the 
line  containing  AY).  Hence  y  can  no  longer  play  the  role  of  the  third  point  in  any  further  solution 
triple  of  the  form  (AZ,Y),  and  we  can  therefore  delete  Y  from  the  end  of  the  active  list.  We  then 
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take  the  dement  Y'  which  becomes  the  new  final  element  of  the  active  list  and  repeat  the  same 
procedure  for  J" .  We  continue  in  this  maimer  until  the  first  time  at  which  the  last  dement  W  on 
the  active  list  is  such  that  the  angle  AWX  is  obtuse.  When  this  happens,  the  triple  (A,X,W)  is  a 
solution  to  Problem  1 ,  because  all  intermediate  points  Q  between  X  and  W  lie  outside  the  cirde 
with  diameter  AX  (since  for  all  these  Q  the  angle  AQX  is  acute).  We  consequently  output  this 
triple,  add  X  to  the  end  of  the  active  list,  and  continue  in  this  manner  until  all  jxjints  in  5  are 
processed. 

It  is  also  easy  to  check  that  the  algorithm  finds  all  solution  triples  to  Problem  1.  Indeed,  let 
(A^,Y)  be  such  a  solution  triple.  It  is  dear  from  the  above  description  that  as  we  process  the  list 
of  points  in  5  sorted  by  their  orientations  about  A ,  every  point  Y  will  be  added  at  some  step  to  the 
top  of  the  active  list.  But  then  the  arguments  in  the  preceding  paragraph  imply  that  Y  will  not  be 
removed  from  this  list  before  X  is  encountered,  and  that  the  processing  of  X  will  eventually  cause 
the  removal  of  each  point  Q  added  to  the  active  list  after  Y  (none  of  which  can  form  an  obtuse 
angle  AQX),  so  that  the  triple  (A^,Y)  will  be  found  by  the  algorithm. 

Finally,  it  is  dear  that  (aside  from  the  sorting  procedure  whidi  is  required  initially)  the 
algorithm  that  we  have  just  described  runs  in  linear  time.  Indeed,  each  of  its  steps  dther  advances 
by  one  item  through  the  sorted  list  of  points  constituting  5,  or  backs  up  along  the  active  list, 
removing  one  element  of  5  at  a  time  from  that  list.  Thus  at  most  2n  steps  will  be  reqiiircd.  Hence 
the  algorithm  is  dominated  by  its  initial  sorting  step,  whidi  consumes  time  0{n  log  n)  for  each 
point  A  (.  S,  and  so  it  runs  in  0(/i'  log  n)  time  overall.  This  condudes  the  proof  of  Lemma  3.2. 
Q.E.D. 

Proof  of  Lemma  3.3: 

The  algorithm  whose  existence  is  asserted  acts  as  foUows  on  each  A  i  S: 

Fix  some  orientation  9  of  the  coordinate  system  (without  loss  of  generality  we  will  refer  to  this  as 
6  =  0,  i.e.  suppose  that  it  is  horizontal.)  Our  basic  idea  is  to  construct  a  list  L  ==  L^  (maintained 
as  a  2-3  tree)  of  all  points  C  in  5  whidi  satisfy  the  following  two  conditions: 
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(1)  C  lies  northeast  of  A,  i.e.,  assuming  A  lies  at  the  origin,  x^  a  0,  y^  ^  O* 

(2)  There  does  not  exist  W  iS  such  that  0  <  y^  <  ^c  and  LAWC  >  y . 

The  algorithm  will  maintain  L  dynamically  while  rotating  the  angle  6  determining  the  orientation 
of  the  coordinate  axes  x,  y  clockwise  through  a  full  360  degrees. 

Note  that  condition  (2)  implies  that  the  part  of  the  interior  of  the  circle  with  diameter  AC 
that  is  cut  off  by  the  horizontal  lines  at  levels  0  and  y^  contains  no  point  W  oi  S.  Note  further 
that  these  two  conditions  imply  that  C  is  minimal  among  all  points  of  5  lying  northeast  of  A, 
under  the  partial  order  (a:  j-)  <  {x'  y)  iff  x  ^  x'  and  y  ^y' .  The  points  in  L  will  be  maintained 
in  ascending  sorted  order  of  their  x-coordinates,  which  means  that  they  are  at  the  same  time 
sorted  in  descending  order  of  their  y  coordinates,  since  for  no  two  clanents  of  such  a  set  of 
minima  do  their  x  coordinates  stand  in  the  same  relationship  as  their  y  coordinates. 


B, 


Bj. 


B,* 


^' 


Hg.  3.6.  An  example.  L  =  [8^,82,8;^^^^^] 
To  see  the  relevance  to  Problem  2  of  the  set  Z.^  of  points  satisfying  conditions  (1)  and  (2), 
suppose  that  for  some  orientation  9  of  the  coordinate  axes  some  point  B  ^  S  comes  to  lie  directly 
to  the  right  of  A.  Let  C  C  L  be  the  rightmost  point  in  L  lying  left  of  S.  Then  we  have 

Lemina  3.5:  If  (A,£,C)  is  a  solution  triple  of  Problem  2  then,  when  B  comes  to  lie  directly  to  the 
right  of  A,  C  satisfies  conditions  (1)  and  (2),  and  is  also  the  rightmost  point  in  L  lying  left  of  8. 
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Proof:  It  is  dear  that,  for  the  orientation  9  of  the  cxx»rdinate  frame  which  makes  B  lie  directly  to 
the  right  of  A,  the  point  C  lies  northeast  of  A.  If  there  existed  a  point  W  i  S  such  that 

0  <  >w  <  yc  ^"'^  LAWC  >  ^,  then  W  would  have  to  lie  in  the  intersection  of  the  horizontal  strip 

{(x  j)  :0  <  y  <  y^}  with  the  open  disc  of  diameter  AC.  But  then  W  would  belong  to  the  interior 
K  of  the  union  of  the  rectangle  ABYX  and  the  two  half-discs  of  diameters  AX,  BY,  as  shown  in 
Fig.  3.3,  contradicting  the  fact  that  (A,B,C)  is  a  solution  of  Problem  2.  Hence  C  satisfies 
conditions  (1)  and  (2).  Moreover,  no  point  D  €  5  lying  right  of  C  and  left  of  B  can  satisfy 
conditions  (1)  and  (2),  for  then  D,  being  minimal,  would  have  to  lie  below  C  and  above  B  (which 
is  at  height  yg  =  0),  again  contradicting  the  fact  that  (A,B,C)  is  a  solution  of  Problem  2.  Q.E.D. 

Hence,  if  as  above  C  C  Z,  is  the  rightmost  point  in  L  lying  left  of  B ,  then  either  the  triple 
(A  fi,Q)  is  a  solution  of  Problem  2,  or  else  no  point  Z?  €  5  can  form  such  a  solution  triple  with 
A,B.  Thus,  if  we  can  maintain  the  sorted  list  L^  efficiently  while  the  coordinate  system  rotates 
clockwise  about  A,  then  the  procedure  just  explained  will  efficiently  find  all  solution  triples 
(A,B,C)  of  Problem  2  (plus  some  spurious  triples). 

Efficient  methods  for  calculating  the  initial  values  of  the  lists  L^,  for  each  A  C  5,  will  be 
described  below.  To  see  what  operations  are  required  for  dynamic  maintenance  of  L^,  note  first 
that  a  new  point  C  ^  S  must  be  inserted  into  Z,^  when,  and  only  when,  the  following  sitiiation 
occurs: 

(i)  C  comes  to  lie  directly  to  the  right  of  A . 

Indeed,  it  is  plain  that,  immediately  after  an  orientation  at  which  this  situation  occurs,  C  will 
belong  to  L^.  Moreover,  if,  for  some  orientation  Oq  of  the  coordinate  axes,  C  lies  northeast  of  A 
but  does  not  satisfy  condition  (2),  it  will  not  satisfy  condition  (2)  for  any  6  lying  clockwise  of  Og 
for  which  C  still  lies  northeast  of  A .   Indeed,  there  exists  W  i  S  such  that  at  Sq  0  <  y^,  <  y^  and 

LAWC  >  —.  Rotating  the  coordinate  axes  clockwise  from  Sq,  these  two  conditions  will  continue 

to  prevail  until  W  comes  to  lie  either  directly  to  the  left  of  A  or  directly  to  the  right  of  C.  But  it  is 
easy  to  sec  that  when  either  of  these  two  cases  occurs,  C  can  no  longer  lie  northeast  of  A.  Hence 
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inserdon  of  C  into  L^  can  only  occur  at  the  orientation  for  which  C  lies  directly  right  of  A,  as 
claimed. 

Next  we  note  that  a  point  C  must  be  deleted  from  L^  only  when  one  of  conditions  (1),  (2) 
comes  to  fail,  i.e.  whenever  one  of  the  following  cases  arises:  either 

(ii)  C  comes  to  lie  directly  above  A,  or 

(iii)  Another  point  W  ^  S,  lying  right  of  A,  comes  to  lie  directly  to  the  left  of  C,  or 

(iv)  Another  point  W  ^  S,  lying  left  of  C,  comes  to  lie  directly  to  the  right  of  A. 

Regarding  case  (iii),  note  that  since  C  a  assumed  to  belong  to  L^  when  such  a  configuration 
arises,  then  the  portion  K  of  the  disc  of  diameter  AC,  bounded  by  AC  and  the  line  passing 
through  C,  W,  contains  no  point  of  S.  Hence  deletions  required  in  case  (iii)  always  correspond  to 
solution  triples  (AJ5,C)  of  Problem  1.  Thus,  by  first  solving  Problem  1,  we  can  prepare  in 
advance  a  superset  Ca>  for  each  A  ^  5,  of  the  set  of  all  pairs  (C,W)  for  which  a  deletion  of  type 
(iii)  from  L^  is  required.  Moreover,  by  Lemma  3.2,  the  total  size  of  all  these  sets  is  at  most 
0(n2  log  n). 

In  regard  to  case  (iv),  we  note  that  it  always  arises  simultaneously  with  the  insertion  of  the 
point  C  into  L,  as  in  (i)  above,  so  that  we  can  combine  deletions  of  type  (iv)  with  corresponding 
insertions  of  tupe  (i)  into  L^. 

To  summarize,  Z-^  is  updated  precisely  as  follows,  as  the  coordinate  system  turns  in  the 
clockwise  direction: 

I.  Whenever  a  point  C  €  S  comes  to  lie  directly  to  the  right  of  A,  insert  C  into  L^  and  ranove  all 
the  points  lying  right  of  C  from  L^. 

n.  Whenever  a  point  C  €  5  comes  to  lie  directly  above  A,  delete  C  from  L^. 

in.  Whenever  two  points  B,C  of  S  become  co-horizontal,  with  B,  C  lying  northeast  of  A,  C  lying 
left  of  B,  and  (A,B,C)  forming  a  solution  of  Problem  1,  delete  B  from  L^  (if  it  is  a  member  of 
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It  remains  to  describe  efficient  techniques  for  calculating  the  initial  values  of  all  the  lists  Z,^, 
for  A  C  5.  For  this  we  again  use  the  solutions  of  Problem  1,  assumed  to  be  available,  and  make 
the  following  observations: 

Let  A,B  d  S,  and  suppose  that  B  d  L^ior  the  initial  orientation  9  =  0  of  the  coordinate  axes,  so 
that  in  particular  B  lies  northeast  of  A.  Let  C  i  S  (resp.  D  e  5)  be  a  point  in  the  upper  half 
(resp.  lower  half)  of  the  disc  of  diameter  AB ,  such  that  the  portion  of  the  disc  bounded  between 
the  diameter  AB  and  the  line  passing  through  B ,  C  (resp.  throxigh  A ,  D)  contains  no  point  of  S  in 
its  interior.  K  such  a  C  exists  and  B  6  Z,^,  then  C  must  plainly  lie  above  B,  and  similarly,  if  such  a 
D  exists  and  B  d  L^,  then  D  must  lie  below  A.  It  is  also  easy  to  sec  that  the  converse  also  holds, 
i.e.  if  B  €  5  is  a  point  lying  northeast  of  A  (for  the  orientation  9  =  0  of  the  axes)  and  if  the 
points  C,  D  respectively  lie  above  B  and  below  A,  then  B  i  L^. 

Thus  for  each  pair  of  points  A,B  ^  S,  with  B  lying  northeast  of  A  when  the  coordinate  axes 
are  given  the  orientation  9  =  0,  we  can  determine  whether  B  i  Lj^  simply  by  examining  the  two 
solution  triples  of  Problem  1  of  the  form  {AJBJ))  and  {B,A,C),  with  D  (resp.  C)  lying  clockwise 
(resp.  counterclockwise)  from  B  (since  we  assume  that  no  three  points  of  5  are  collinear,  these 
solution  triples  are  unique  if  they  exist),  and  checking  that  C  lies  above  B  and  that  D  lies  below  A. 

It  is  plain  from  Lemma  3.2  that  the  total  time  required  by  the  procedure  that  we  have 
outlined  is  0{rr  log  n). 

Summarizing,  we  see  that  the  following  algorithm  efficiently  calaJatcs  a  superset  of  the  set 
of  solution  triples  to  Problem  2: 

(1)  Using  the  solution  to  Problem  1,  construct  (in  time  C>(n^og  «))  for  each  A  €  J  the  list  C^, 
of  coordinate-axis  orientations  9  for  which  there  exist  B,C  i  S  which  are  co-horizontal  at  9, 
and  which  with  A  form  a  solution  triple  (C,A  ,5)  of  problem  1. 

(2)  Again  using  the  solutions  to  Problem  1,  employ  the  procedure  described  above  to  compute 
the  initial  values  of  all  the  lists  L^  for  the  orientation  9  =  0  of  the  coordinate  axes  and  for 
all  A  €  5.  This  requires  time  0(n-  log  n). 
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(3)  For  each  A  €  5  sort  the  2n  +  \Q^\  critical  orientatiom  0se<2iT  which  either  belong  to  Q^, 
or  are  orientations  at  which  some  B  ^  S  comes  to  lie  directly  to  the  right  of  A  or  directly 
above  A .  These  orientations  should  be  sorted  into  increasing  counterclockwise  order. 

(4)  Process  the  critical  orientations  produced  by  step  (3)  in  counterclockwise  order.  For  each 
such  orientation  9  do  the  following. 

(4. a)  If  at  e  some  5  C  5  lies  directly  to  the  right  of  A ,  find  the  rightmost  point  C  in  L  lying  left  to 
B,  and  output  the  triple  (A^,C)  if  |A8|  ^  d.  In  addition,  insert  B  into  L  and  delete  all 
points  lying  right  to  B  from  L.  (This  and  the  two  substq»  described  just  below  can  be 
accomplished  in  0(log  n)  time,  if  L  is  maintained  as  a  2-3  tree.) 

(4.b)  If  at  8  some  fl  €  5  lies  direcdy  above  A,  delete  B  from  (the  beginning  of)  L. 

(4.c)  If  9  €  Oxt  i-c.  if  two  points  B,C  ^  S  become  co-horizontal  at  9  and  (A,C,B)  solves 
problem  1,  delete  C  from  L  if  C  6  L;  otherwise  do  nothing. 

The  correctness  of  the  algorithm,  in  the  sense  that  it  produces  all  correct  solution  triples, 
plus  some  spurious  triples,  follows  from  the  preceding  arguments.  The  time  required  by  the 
algorithm  is 


0(  2  (2«  +  ICaI)  log  n  )  =  Oin^og  n) 
AiS 


as  asserted.  Q.E.D. 


Remark:  It  is  worth  noting  one  last  remark  concerning  the  "harmlessness"  of  the  spurious  triples 
produced  by  the  algorithm.  Suppose  that  (AJB,C)  is  such  a  spurious  triple.  Then  the  "racetrack" 
shown  below  contains  a  point  £>  6  5  in  its  interior. 
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Fig.  3.7.  A  spurious  triple  (A,B,C) 
Suppose  that  we  apply  the  previously  described  procedure  for  calculating  the  associated  Voronoi 

vertices  of  types  Il.a  and  n.b  to  the  pair  A,B  d  comers.  Because  of  the  presence  of  B,  the 
interval  L(y)  of  placements  available  for  the  ladder  at  height  y  above  AB,  will  shrink  to  a  jxjint 
and  vanish  at  a  height  y^  smaller  than  yQ  =  half  the  distance  of  C  from  AB .  Hence  this  procedure 
will  terminate  before  reaching  the  critical  height  y^  associated  with  C,  and  therefore  will  not  be 
affected  by  the  presence  of  this  spurious  critical  height. 

Next  we  proceed  to  complete  the  analysis  of  Problem  3,  i.e.  to  prove  Lemma  3.4. 

Proof  of  Lemma  3.4:  For  each  point  A  ^  S  let  B  d  S  be  another  point  whose  distance  from  A  b 
^d.  For  each  such  B  let  B^  denote  the  point  at  distance  d  from  A  which  lies  on  the  ray  from  A 
through  B,  and  let  5j(A)  denote  the  collection  of  all  these  "displaced"  points  B^.  Suppose  that 
(A,B,C)  is  a  solution  triple  to  Problem  3,  let  AB^YZ  be  the  corresponding^ectangle,  and  let  H'  be 
the  midpoint  of  B^Y.  Then  W  must  lie  on  the  standard  Voronoi  diagram  Voro(5U5j(A))  of  the  set 
S\JS^(A).  Indeed,  the  object  in  J  nearest  to  W  is  C,  by  definition,  and  W  is  also  at  the  same 
distance  r  from  Bj.  However,  no  other  point  K^  in  S^(A)  is  nearer  to  W  than  B^,  for  if  it  were,  K 
would  have  to  lie  both  inside  the  circle  of  radius  r  about  W  and  the  drdc  of  radius  d  about  A ; 
hence  inside  the  rectangle  AB^YZ  which  is  impossible.  Thus  W  must  lie  on  the  boundary  of  the 
Voronoi  cell  V(Bj)  consisting  of  all  points  nearer  to  B^  than  to  any  other  object  in  5U5j(A).  This 
cell  is  easily  seen  to  be  star-shaped  with  respect  to  B^;  hence  any  line  through  B^  intersects  its 
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boundary  in  just  two  points.  These  observations  justify  the  following  algorithm: 

For  each  comer  A  in  S,  build  the  set  S^(A)  and  compute  the  standard  Voronoi  diagram  of 
SUS^(A).  Then,  for  each  B^  €  S^iA)  draw  the  line  through  B^  perpendicular  to  the  segment  AB^, 
and  find  the  two  intersections  of  this  line  with  the  boimdary  of  the  cell  V(B^)  in  VorQ(SU5^(A)). 
These  intersections  define  all  the  rectangles  which  might  correspond  to  solutions  of  Problem  3. 
Since  the  standard  Voronoi  diagram  of  SUS^(A)  can  be  constructed  in  time  0(n  log  n)  (using  a 
recent  algorithm  by  Yap  [Ya];  see  also  [Ki]),  while  the  search  for  intersections  of  this  Voronoi 
diagram  with  the  perpendicular  lines  just  mentioned  can  be  performed  in  time  proportional  to  the 
number  of  edges  of  the  diagram,  i.e.  in  0(n)  time,  it  follows  that  the  time  required  to  process 
each  comer  A  ^  S  is  0(n  log  n).  Repeating  this  procedxire  for  each  A  ^  S  gives  an  0(n^  log  n) 
algorithm  which  generates  a  set  of  size  0(n^  inrluding  all  solution  triples  to  our  problem.   Q.E.D. 

4.  Voronoi  Verticea  of  Type  n.c 

Next  we  estimate  the  number  of  Voronoi  vertices  of  type  n.c  and  describe  an  efficient 
algorithm  for  their  calculation. 

Theorem  4.1:  There  are  at  most  0(n^  log'n)  Voronoi  vertices  of  type  n.c,  and  they  can  all  be 
calculated  in  time  0(n^og  n  log'n). 

Proof:  Let  s^  and  ^2  be  two  comers  in  5,  and  consider  the  set  of  Voronoi  vertices  of  type  n.c 
associated  with  s^  and  ^2.  i-c-  placements  Z  in  which  s^  (resp.  s^)  touches  the  left  (resp.  the  right) 
side  of  the  (clear)  R(Z),  while  two  other  objects  toiich  opposite  bends. 

Let  D  be  the  midpoint  of  the  segment  s^S2.  Rotating  and  translating  axes  if  necessary,  we  can 
assume  that  £>  is  the  origin  and  that  j^,  13  Uc  on  the  x-axis  at  distance  2a  apart,  with  s^  to  the  left 
of  j2.  For  any  Voronoi  vertex  Z  of  type  H.c  associated  with  j^  and  32,  D  must  lie  on  B(Z)  and  Z 
must  have  an  orientation  which  lies  strictly  between  0  and  ti.  Since  D  is  on  5(Z),  Z  must  have  the 
form  (r  cos  9,  r  sin  9,  9)  for  some  r  and  9,  and  the  distances  from  B(Z)  to  s^  and  ^2  niust  both  be 
a  sin  9. 
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For  each  orientation  0<6<ir  of  B  let  L(9)  denote  the  set  of  all  r<0  such  that 
R(Z,  a  sin  9),  where  Z  =  (r  OM  9,  r  sin  9,  9),  is  dear  and  touches  s^  on  its  left  side  and  j,  on  its 
right  side.  If  the  set  L(9)  is  nonempty  then  it  is  easily  seen  that  it  is  an  interval  [/■(9),  g(9)],  where 
/(9)  and  ^(9)  are  both  negative.  In  this  case  wc  can  let  a(9)  (resp.  p(9))  denote  the  set  of  objects 
of  S  meeting  the  lower  (resp.  upper)  bend  of  R(f(Q))  (resp.  /?Cs(9))). 

In  contrast  to  the  preceding  analysis,  in  which  i4(f)  changed  monotonically  with  f,  the  interval  L(Q) 
need  not  be  monotone  in  9.  Nevertheless  we  can  track  the  way  in  which  it  varies  in  much  the 
same  way  as  before.  Specifically,  we  note  that,  as  we  vary  the  orientation  9,  the  function  g(Q) 
remains  continuous  and  the  set  p(9)  stays  constant,  except  at  "critical"  orientations  9  (sec  Fig.  4.1) 
at  which  one  of  the  following  conditions  holds: 
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Fig.  4.1.  Critical  orientations  in  the  calailation  of  type  n.c  vertices 
(c.l)       Two  objects  s^,  Si^S  meet  the  upper  bend  of  the  racetrack  R  =  R(s(9)).   (In  this  case 

the  set  P(e')  can  change  as  9'  passes  through  9.) 
(c.2)  A  third  comer  Jj  €  5  touches  one  of  the  sides  of  R,  so  that  the  projection  of  j,  onto  the 
ladder  lies  above  the  projections  of  the  other  two  comers  s^  and  J2-  In  addition,  dthcr  a 
fourth  object  J4  touches  the  upper  bend  of  /?  (in  which  case  contact  with  s^  defines  the 
positional  limit  ^(9)),  or  one  of  the  comers  jj,  53  projects  onto  P  (in  which  case  this 
relationship  defines  g(Q).  (In  this  case  the  ladder  must  be  moved  down  to  the  left  until  J3 
contacts  its  upper  bend  in  order  to  continue  to  turn  counterclockwise  thitnigh  the  critical 
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orientation  9.  Thus  g(9')  change*  discontinuoiisly  at  9'  =  9.) 
(c.3)  A  third  cxjmer  jj  €  S  touches  one  of  the  sides  of  R  and  the  projection  of  53  onto  the 
ladder  lies  between  the  projections  of  the  two  comers  jj  and  J2-  1°  addition,  either  a 
fourth  object  s^  touches  the  upper  bend  of  R  (in  which  case  contact  with  s^  defines  ^(9), 
or  one  of  the  comers  s^,  S2  projects  onto  P  (in  which  case  this  relationship  defines  ^(9). 
(In  this  case  the  ladder  cannot  continue  to  turn  counterclockwise  and  remain  nearest  to  j^ 
aiul  J2  while  passing  through  D,  so  that  Z.(9')  becomes  empty  for  9'  slightly  greater  than 
9.) 

This  last  condition  implies  the  following  somewhat  more  convenient  condition,  which  is  better 
suited  for  calculation  of  the  critical  orientations  6: 

(c.3')  The  comer  J3  lies  in  the  upper  semicircle  C  whose  diameter  is  j^jji  ^^  ^  section  of  C 
bounded  by  s^S2  and  the  ray  s^s^  (or  by  3^32  and  the  ray  jj^s)  >*  dear.  (Indeed,  in  an 
extreme  configuration,  J3  would  line  up  with  one  of  3^  and  jj  (say  32  for  deflniteness)  to 
form  a  right  triangle  3^3^32  (see  Fig.  4.1  (C.3)),  with  J3  exactly  on  the  scmidrde  whose 
diameter  is  3^32.  In  all  other  configurations  of  type  (C.3),  53  lies  within  this  semidrde; 
and  the  pxjrtion  of  the  semidrde  cut  off  by  the  line  3^3y  (which  line  is  a  side  of  the 
racetrack)  is  easily  seen  to  be  dear.) 

Finally,  the  placement  Z  =  (g(9)  cos  9,  g(9)  sin  9,  9)  may  be  a  Voronoi  vertex  of  type  n.c  that 
we  seek.  (In  this  case/(9)  =  g(9),) 

Similar  and  symmetric  conditions  govern  the  oocxirrence  of  discontinuous  changes  in  /(9)  and 

Qt(9). 

Let  N(3^,3^  be  the  number  of  critical  orientations  satisfying  one  of  the  conditions  (c.l)-(c.3') 
(or  one  of  their  symmetric  counterparts),  and  suppose  that  aD  these  critical  orientations  have 
already  been  calculated.  Then  in  each  of  the  N{3^^2)  "•■  1  intervals  of  orientations  bounded  by  two 
adjacent  critical  orientatiom  (or  by  a  single  critical  orientation  and  by  one  of  the  eoc^>oina  0  or  it 
of  the  available  range  of  orientations)  the  two  functions  /(9)  and  g(9)  vary  continuously  and  the 
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(singleton)  sets  a(e)  and  p(9)  rranain  constant.  Given  a(9)  and  P(8)  we  can  express /(e)  and  g(6) 
by  entirely  elementary  formulae  and  hence  determine  in  constant  time  whether  such  an  interval  / 
contains  orientations  6  for  which  /(9)  =  g(9)  corresponding  to  some  Voronoi  vertex  of  type  n.c, 
i.e.  to  a  critical  orientation  at  which  f  =  g.  Note  that  in  such  an  interval  there  can  be  at  most 
0(1)  9  satisfying /(9)  =  g(9)  since  both/  and  g  will  be  elanentary  analytic  functions  in  /.)  Thus, 
once  the  critical  orientations  are  available,  we  can  find  the  type  n.c  vertices  (associated  with  9 
satisfying  /(9)  =  ^(9))  in  time  0(Nis^,s^  log  N(si^^). 

Hence  it  only  remains  to  show  how  the  sets  of  critical  9  satisfying  one  of  the  conditions  (c.l) 
-  (c.3')  can  be  found  in  time  Oin^og'n).  For  this,  we  argue  as  follows.  Orientations  9  satisfying 
condition  (c.l)  correspond  to  Voronoi  vertices  of  type  I,  of  which  there  are  at  most  0{n^  log*n), 
and  which  have  all  been  found  already.  For  each  placement  at  which  condition  (c.2)  (resp.  (c.3')) 
holds,  one  of  the  triples  (3^,3^^^  or  (^2,^3^1)  is  clearly  a  solution  to  Problem  2  (resp.  Problem  1) 
of  Section  3  above.  Hence  in  total  (i.e.  summed  over  all  pairs  of  comers  s^,  $2  i  S)  there  are  at 
most  0(n^  critical  orientations  of  either  of  types  (c.2)  or  (c3'),  and  they  can  easily  be  calculated 
from  the  solution  to  Problems  1  and  2.  (Recall  that  the  algorithm  given  in  Section  3  for  solving 
Problem  2  produces  a  superset  of  the  set  of  all  true  solutions  to  that  problem.  However,  arguing 
as  in  Section  3,  we  can  see  that  the  spurious  triples  produced  by  the  algorithm  will  not  cause  any 
problem,  because  they  can  easily  be  detected  and  discarded  as  they  are  encountered  during  the 
progressive  inaease  in  9.) 

All  this  implies  that  there  are  altogether  at  most  0{n^  log'n)  critical  orientations  of  types 
(c.l)-(c.3'),  summed  over  all  pairs  of  comers  j^,  ^2  ^  S,  and  thus  also  at  most  0{n^\og*n) 
Voronoi  vertices  of  type  n.c.  Thus  altogether  the  total  time  required  to  calculate  these  vertices  is 

2  ^(*i.*2)  log//(5i,J2)  =  0(n2log  n  log'n) 
as  asserted.  Q.E.D. 

Remark:  In  the  preceding  descriptions  of  the  various  procedures  for  the  calculation  of  free 
Voronoi  vertices  we  have  tadtly  assumed  that  the  four  objects  s^,  $2,  s^  and  s^  defining  such  a 


-38- 

vertex  arc  separated.  In  fact,  as  noted  in  the  preceding  paper  of  this  series  (see  [OSY]),  some  of 
s^,  S2,  53,  Si  can  be  adjacent.  Nevertheless  the  procedures  that  we  have  developed  can  easily  be 
adapted  to  handle  nonseparated  objects.  To  illustrate  the  issues  involved,  consider  Voronoi 
vertices  of  type  I  for  which  the  end  P  of  the  ladder  is  simultaneously  nearest  to  two  nonseparated 
objects,  e.g.  a  wall  edge  s^  and  one  of  its  comers  J2-  Then  P  is  easily  seen  to  lie  on  the  line 
perpendicular  to  s^  at  J2-  ^  >°  f^li  ^^  ^"^  <=^>Q  ^  regarded  as  part  of  the  standard  Voronoi 
diagram  VorQ^S),  making  it  possible  to  use  the  analysis  of  Section  2  verbatim  in  this  case  too. 

Next  consider  vertices  Z  of  type  n.a  or  n.b  in  which  a  wall  edge  j^  and  one  of  its  comers  Jt  lie 
parallel  to  B(Z)  with  J2  nearest  to  an  interior  point  of  B(Z)  (this  corresixmds  to  the  case 
considered  in  Section  3,  where  s^  and  ^2  are  assumed  to  be  two  separated  comers).  The  present 
case  is  even  simpler  than  that  treated  in  Section  3,  because  at  most  0(n)  pairs  s^,  Sj  need  be 
considered,  making  it  much  easier  than  in  Section  3  to  compute  and  process  the  relevant  critical 
heights.  We  leave  it  to  the  reader  to  verify  that  the  associated  vertices  of  these  types  can  be  found 
in  time  0(n^  log  n).  Finally,  we  observe  that  verticra  of  type  II.c  occur  only  when  the  ladder  is 
nearest  to  four  separated  objects,  and  no  cases  of  this  type  involving  non-separated  objects  arise. 

The  following  theorem  summarizes  the  preceding  discussion: 
Theorem  4.2:  The  total  number  of  Voronoi  vertices  in  Vor^(S)  is  0(n-  log'n)  and  they  can  all  be 
calculated  in  time  0(n^  log  n  log'n). 

5.   Computation  of  Semi-free  and  Pseado- Vertices. 

Now  that  we  have  calculated  all  the  "true"  Voronoi  vertices,  we  can  go  on  to  consider  the 
problem  of  calculating  all  the  auxiliary  semi-free  and  pseudo  vertices  whidi  appear  in  the  skeleton 
N  in  which  we  are  ultimately  interested. 

Elfidcnt  calculation  of  these  auxiliary  vertices  is  much  simpler  than  calculation  of  true 
vertices  as  described  in  the  preceding  sections,  because  each  pseudo  and  semi-free  vertex  is 
determined  by  three  rather  than  four  objects  of  5.  Moreover,  most  of  the  methods  presented  in 
the  preceding  sections  for  the  calculation  of  true  Voronoi  vertices  can  be  adapted  to  calculation  of 
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scmi-free  and  pseudo  vertices. 

We  begin  by  considering  semi-free  Voronoi  vertices.  It  is  easily  checked  that  only  four  types 
of  semi-free  vertices  need  be  considered: 

(A)  Vertices  Z  at  which  P  touches  a  comer  j  €  5  (cf .  Fig.  5.1(a)). 


(a) 


'^^^777^777^^777^777777- 


(b) 


(c) 


Q 


(d) 


Fig.  5.1 

(B)  Vertices  Z  for  which  B(Z)  touches  one  edge  of  5  at  P,  another  edge  at  Q,  and  some  comer  of 
5  at  an  interior  point  on  B  (cf.  Fig.  5.1(b)). 

(Q  Vertices  Z  for  which  fl(Z)  touches  one  edge  of  S  at  P,  and  two  comers  of  S  at  two  interior 
points  on  B,  where  these  two  comers  of  5  lie  on  the  same  side  of  B  (cf.  Rg.  5.1(c)). 

(D)  Same  as  (Q,  except  that  the  two  comers  of  S  that  B  touches  lie  on  different  sides  of  B. 

Type  (A)  semi-free  vertices  are  easy  to  find:  Put  P  at  each  of  the  0(n)  comers  of  5,  and  for 
each  such  comer  s  find  all  the  orientations  9  for  which  fl(j,9)  touches  another  object  of  5,  but 
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does  not  penetrate  any  of  these  objects.  This  can  easily  be  done  in  0{n  log  n)  time,  using  a  simple 
sweeping  technique. 

To  find  vertices  of  types  (B)-(D),  we  can  use  a  procedure  quite  similar  to  that  used  for  true 
Voronoi  vertices  of  type  I.  However,  instead  of  moving  P  along  a  Voronoi  edge  of  Vor(j(5),  wc 
move  it  along  one  of  the  0{n)  edges  s  d  S,  while  maintaining  a  list  L  of  the  critical  orientations  9 
at  v^di  B  touches,  but  does  not  penetrate  another  object  of  S.  Straightforward  analysis  of  the 
situation  shows  that  L  needs  to  be  updated  only  at  0{n)  positions  along  s,  and  that  the  semi-free 
vertices  which  we  seek  correspond  to  positions  at  which  two  such  critical  orientations  coincide,  and 
so  can  be  obtained  as  L  is  being  updated.  Although  some  of  the  pertinent  details  are  not 
immediate,  none  are  difficult,  and  we  therefore  leave  their  elaboration  to  the  reader. 

With  these  provisos  we  obtain 

Theorem  5.1:  There  are  at  most  0{n^  semi-free  Voronoi  vertices,  all  of  which  can  be  found  in 
time  0(/i^  log  n). 

Next  we  turn  to  the  computation  of  pseiido  vertices.  Recall  that  these  are  placements 
Z  =  (x^y,9)  of  B  such  that  fl(Z)  is  nearest  to  three  objects  Ji,  S2,  and  jj  of  5,  with  j^  and  ^2 
separated,  and  such  that  Z  lies  on  the  core  streamline  of  the  midsurfaoe  between  j,  and  s^  (sec 
[OSY]);  that  is,  the  orientation  9  of  Z  is  horizontal  relative  to  the  canonical  position  of  jj  and  Ji 
(see  the  paragraph  preceding  Fig.  2.1  for  a  description  of  these  canonical  positions).  As  in  Section 
2,  we  shall  refer  to  the  three  possible  combinations  of  types  of  jj  and  Ji  ^  'wall-wall',  Vall- 
comer',  and  'comer-comer'. 

Rrst  consider  the  wall-wall  case  (see  Fig.  2.1).  It  is  easily  seen  that  the  Voronoi  pseudo 
vertices  Z  in  qxiestion  are  such  that  P  lies  on  the  standard  Voronoi  diagram  Voro(S).  Hcnoe  there 
are  at  most  0{n)  possible  pairs  {s'^^-i)  of  walls  that  need  be  considered.  Let  {s^^yS^  be  one  of  these 
pairs,  and  let  «  be  an  edge  of  Voro(S)  consisting  of  points  simultaneously  dosest  to  j.  and  s-,-,  this 
edge  plainly  lies  on  the  bisector  L  of  the  angle  between  s^  and  i-,.  For  each  additional  object  s  i  S 
we  can  find  in  constant  time  all  horizontal  placements  Z  of  B  at  which  /*  lies  on  *  and  fl(Z)  is 
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equidistant  from  jj,  S2,  and  s.  Sorting  these  0(n)  placements  by  the  positions  of  P  along  e,  and 
processing  them  in  increasing  order,  it  is  easy  to  find  all  placements  which  are  pseudo  vertices  of 
the  type  we  seek,  in  time  0(n  log  n).  Since  there  are  0(n)  edges  in  Voro(5),  it  follows  that  at 
most  0(n^  pseudo  vertices  can  arise  in  the  wall-wall  case,  and  that  they  can  all  be  found  in  time 
0(n^  log  «). 

The  comer-comer  case  can  be  handled  in  exactly  the  same  manner,  due  to  the  fact  that  in 
this  case  each  corresponding  pseudo  vertex  {X,0)  must  once  again  represent  a  ladder  position  lying 
along  the  line  exactly  half-way  between  the  vertically  separated  comers  s^  and  S2. 

The  wall-comer  case  requires  a  slightly  modified  treatment.  Suppose  that  *i  is  a  comer  and 
J2  a  wall  edge  in  their  canonical  position  (i.e.  J2  00  the  x-axis,  s^  on  the  positive  ^-axis),  and  let 
Z  =  (X,0)  be  a  pseudo  vertex  nearest  to  s^,  S2  and  some  third  object  sy  Then  there  are  two 
possible  subcases:  either  one  of  the  endpoints  i"  or  C  of  the  ladder  lies  on  Voro(S),  in  which  case 
much  the  same  procedure  as  has  just  been  sketched  can  be  applied,  or  B(Z)  a  parallel  to  ^2  and 
passes  throtigh  the  midpoint  X  of  the  segment  s^A,  where  A  is  the  foot  of  the  perpendicular  from 
jj  to  ^2-  In  this  latter  case  the  point  X  must  lie  on  VorQ(S),  and  again  much  the  same  procedure  as 
has  just  been  sketched  can  be  called. 
The  following  Theorem  summarizes  the  preceding  discussion: 

Theorem  5.2:  There  are  0(n^  Voronoi  pseudo  vertices,  and  they  can  all  be  calculated  in  time 
Oin^  log  n). 

This  concludes  our  description  of  the  procedures  for  the  calculation  of  the  vertices  of  the 
Voronoi  skeleton  N.  Our  analysis  of  N  will  be  completed  in  the  next  section,  which  explains  how 
relationships  of  vertex  adjacency  along  Voronoi  edges  and  streamlines  of  the  skeleton  N  can  be 
determined. 

6.  Adjactncy  rdattonshlps  fai  the  skeleton  N. 

In  this  section  we  move  on  to  the  final  task  required  for  construction  of  the  Voronoi  diagram 
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for  a  ladder,  namely  we  find  all  the  Voronoi  edges  and  streamlines  cxmnccting  the  prcviouily 
calculated  Voronoi  vertices,  pseudo  vertices  and  semi-free  vertices.  Since  the  total  number  of  these 
vertices  has  been  shown  to  be  0{rP-  log'n),  and  since  eadi  vertex  of  VorJ^S)  is  an  cndpoint  of  at 
most  four  edges  (sec  [OSY]),  it  follows  easily  that  the  total  number  of  edges  in  the  diagram  is  also 
0{n'  log'n).  To  actually  determine  adjacency  of  vertices  by  finding  these  edges  and  streamlines, 
we  can  use  the  following  straightforward  procedure: 

(1)  For  each  (free,  semi-free  or  pseudo)  vertex  Z  we  can  determine  (in  constant  time)  the  few 
virtual  edges  (called  midlines  in  [OSY])  and  streamlines  on  which  Z  lies.  As  is  shown  in  [OSY], 
each  connected  portion  of  any  of  these  1 -dimensional  curves  can  be  parametrized  in  a  standard 
way  by  a  parameter  t  ranging  over  the  dosed  interval  [0,1]  or  the  unit  cirde  5\  and  for  each 
placement  Z  lying  on  such  a  curve  the  value  of  the  parameter  t  corresponding  to  Z  can  be 
calculated  in  constant  time. 

Thus  (in  overall  time  0(n^  log  n  log'n))  we  can  collect  the  set  of  all  vertices  Z  lying  on  any 
midline  or  streamline  "y,  and  can  sort  these  placements  in  increasing  order  of  the  corresponding 
parameter  f.  Note  however  that  some  placements  Z  lying  on  such  a  "y  may  not  actually  lie  on  the 
Voronoi  diagram  Vor^(S),  because  fl(Z)  might  be  nearest  to  some  object  in  S  other  than  those 
objects  defining  the  midline  or  streamline  -y.  However,  the  points  along  7  at  which  this  first  starts 
to  happen  are  dearly  Voronoi  vertices.  Thus  we  need  to  find  all  the  subarcs  of  y  delimited  by  two 
such  Voronoi  vertices.  Let  the  parametrization  of  -y  be  denoted  by  W(t),  t  €  [0,1],  and  let  K^ 
denote  the  set  of  all  vertices  lying  on  -y  sorted  b  increasing  order  of  r.  Given  any  Z  =  W^t^)  €  K^ 
and  its  corresponding  parametric  value  f^,  we  can  easily  determine  whether  for  t  slightly  near  t^ 
placements  W(t)  €  -y  actually  belong  to  a  Voronoi  edge.  This  can  be  done  by  computing  the 
distance  of  placements  belonging  to  W(t)  to  the  objects  nearest  the  placement  W^t^)  for  each 
infinitesimal  deviation  of  t  from  r^).  Once  this  information  is  available,  all  that  remains  to  be 
done  is  to  find  all  pairs  Zj,  Zj  of  adjacent  vertices  b  K^  (with  r^  <t^  such  that  for  t  between  r^ 

and  I7  the  placements  W{i)  lie  on  a  Voronoi  edge.  Qearly  any  two  such  placements  Zj  and  Z;  are 
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adjacent  to  one  another  in  the  skeleton  N,  and  moreover  these  are  the  only  adjacent  pairs  of 
vertices. 

This  gives  us  the  skeleton  N.  The  time  required  by  this  procedure  to  accomplish  initial 
distribution  of  the  Voronoi  vertices  into  the  various  sets  K^,  and  to  sort  the  vertices  in  each  of  the 
sets  K  in  parametric  order,  is  plainly  0{n^  log  n  log'n). 

Finally,  we  note  that  any  placement  Z  i  FP  can  be  retracted  to  N  using  straightforward 
computations  requiring  much  less  than  0(n^  time  (see  [OSY]  for  more  detail).  Thus  the 
retraction-based  motion  planning  algorithm  for  a  ladder  that  results  from  the  preceding  analysis 
will  run  in  0(n^  log  n  \og*n)  time,  as  daimed. 

Remark:  To  get  some  sense  of  how  far  the  algorithm  that  we  have  sketched  deviates  from  best 
possible  performance,  it  is  worth  noting  that  in  the  worst  case  the  space  FP  of  free  placements  of 
B  can  have  Cl(n^  connected  components.  To  see  this,  consider  the  situation  shown  in  Fig.  6.1,  in 
which  the  free  space  V  is  bounded  by  a  convex  polygon  and  contains  two  parallel  rows  of  n  small 
obstacles  each.  If  the  dimensions  of  the  enclosing  polygon  and  the  locations  of  the  interior 
obstacles  are  chosen  appropriately,  the  ladder  B  can  be  placed  between  any  pair  of  adjacent 
obstacles  in  the  i^)per  row,  and  between  any  pair  of  adjacent  obstacles  in  the  lower  row,  while 
there  is  no  way  to  move  B  from  such  a  placement  to  another  placement  between  two  other  pairs 
of  adjacent  obstacles,  making  it  plain  that  in  this  case  the  space  FP  of  free  placements  has  C^n?) 
connected  components. 
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Fig.  6.1 
Condndlng  ranurk:  All  the  algorithim  given  in  Secdom  3  and  4  run  in  time  0((n^  +  K)  log  n), 
where  AT  is  the  number  of  Voronoi  vertices  of  type  L  Hence  the  factor  log'/i  appears  in  the 
preceding  complexity  estimates  only  because  it  appears  in  our  estimate  of  the  number  of  type  I 
Voronoi  vertices,  as  given  in  Section  2.  It  therefore  seems  quite  plausible  that  a  more  detailed 
analysis  of  the  behavior  of  the  functions  X,(9)  defined  in  Section  2  can  improve  the  bounds 
established  on  the  size  of  Vor^(S)  and  hence  on  the  total  running  time  of  the  algorithms  presented 
in  this  paper. 
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